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foreword 


This quarterly review of reactor development has been prepared at the request 
of the Division of Information Services of the United States Atomic Energy Com- 
mission. Its purpose is to assist interested organizations in the task of keeping 
abreast of new results in reactor technology for civilian application. 

The report is a concise discussion of selected phases of research and develop- 
ment for which there have been significant advances or a heightened interest in 
the past few months. It is not meant to be a comprehensive abstract of all mate- 
rial published during the quarter nor is it meant to be a treatise on any part of 
the subject. The intention is to cover the various areas of reactor development 
from the general viewpoint of the reactor designer rather than from the more 
detailed points of view of specialists in the individual areas. However, papers 
which are thought to be of particular significance or of particular usefulness in 
specialized fields will be mentioned in short notes. In the over-all plan of the re- 
port, it is intended that various subjects will be treated from time to time and 
will be brought up to date at that time. 

Any interpretation of results which is given represents only the opinion of the 
editors of the report, who are General Nuclear Engineering Corporation person- 
nel. Readers are urged to consult the original references in order to obtain all of 
the background ofthe work reported and to obtain the interpretation of the results 
given by the original authors. 


W. H. ZINN 
General Nuclear Engineering Corporation 
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APPLICATIONS: NUCLEAR SHIP PROPULSION 





General Economic Considerations 


The success of nuclear submarines in the United 
States has demonstrated the technical feasibility 
of nuclear ship propulsion. The economic feasi- 
bility in commercial applications remains to be 
demonstrated. A number of reports, one cover- 
ing a very extensive study,' have recently treated 
the subjects of nuclear ship propulsion and 
its economic relationships. 

Many of the economic factors in ship propul- 
sion are similar to those which govern the 
economics of central station nuclear power. 
Two new variables enter, however: (1) the 
weight of the nuclear plant per unit of output 
(which in comparative evaluations is balanced 
against the weight of a conventional propulsion 
plant plus that of its fuel requirement for a 
given voyage), and (2) the speed at which the 
ship is to be driven (which is determined by 
the capacity of the plant installed in the ship). 
In addition to these new considerations, the 
consideration of reactor safety enters the eco- 
nomic picture perhaps more strongly than in 
the case of the central station reactor because 
(1) the propulsion reactor is more subject to 
external hazards, such as collisions, and (2) be- 
cause the hazard problems are international 
rather than national ones. 

Aside from the foregoing considerations, the 
main economic characteristics, as in the case 
of the central station plant, can be grouped 
into those which affect capital cost and those 
which affect fuel cost. There is no doubt that 
the capital cost of a nuclear propulsion system 
will exceed that of a conventional system by a 
very substantial margin. With regard to capital 
cost the nuclear ship propulsion plant has, in 
one respect, a less favorable competitive situa- 
tion than does the central station reactor, for 
there does not appear to be an immediate 


prospect for very high-power propulsion plants, 
and it is well established that the nuclear 
plant competes more favorably in the larger 
sizes. In another respect, however, its position 
is more favorable. Its cost represents only a 
fraction of the total investment in the ship, and 
therefore the economic effects of performance 
improvements, such as higher speed or higher 
use factor, which the nuclear plant may provide, 
are proportionately magnified. 

In the area of fuel costs nuclear plants are 
potentially cheaper than conventional plants, 
and indeed, it may be said that unless this 
potentiality is realized there is little hope of 
competing successfully. The fuel oil used for 
marine fuel is considerably more expensive than 
the fuel used in the average central station 
plant. Reference 1 quotes a value of $3.00 
per barrel as an average price at the docks in 
the United States; this corresponds to about 
50 cents per billion British thermal units or 
2 mills per thermal kilowatt hour. However, 
the attainment of low fuel cost will be more 
difficult in propulsion reactors than in central 
station reactors because of the high neutron 
loss by leakage which tends to characterize 
the relatively small propulsion reactors. Fur- 
ther difficulties may arise in specific cases 
because of the very great emphasis on reactor 
safety. For example, this emphasis may dictate 
the use of steel rather than zirconium for fuel 
jackets in water-cooled reactors. 

Because of the performance variables of ship 
speed and specific weight of the propulsion 
plant, as well as the variety of ship applica- 
tions, there exists a wide range of possible 
circumstances under which nuclear propulsion 
might be considered, and there is reason to 
expect that nuclear propulsion will be competi- 
tive in specific applications long before it be- 
comes generally competitive. 
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Current ship designs and design speeds have 
evolved from the economic relations peculiar to 
their fossil-fueled propulsion plants. It would 
be expected that the application of nuclear 
power would lead to changes inthese character- 
istics which would improve the competitive 
position of the nuclear plant over that which 
would apply if a nuclear plant were simply 
substituted for a conventional plant in a specific 
ship. In particular, it would be expected that 
the speeds of nuclear ships would be higher, 
since the main restrictions on speed in con- 
ventional ships are the cost and the weight and 
space requirements of fuel oil. However, the 
speed-resistance relationship for surface ships 
will probably limit the speed increase to a 
modest one. For a given ship, the required 
propulsion power increases with speed (V) at a 
rate higher than the third power of V. Con- 
sequently, in current studies of nuclear pro- 
pulsion for nonpassenger ships, one does not 
find postulated speeds higher than about 25 
knots. For completely submerged vessels the 
power requirement increases about v*. hence, 
for such types, higher speeds might be con- 
sidered. Reference 2, which treats a submarine 
tanker of 30,000 deadweight ton cargo capacity 
and a cruising speed of 22 knots, finds that at 
speeds less than about 19 knots the submarine 
power requirements are higher than those of a 
surface tanker of equal capacity, but at higher 
speeds the surface propulsion requirements be- 
come much higher. For example, at the cross- 
over point of 19 knots the required power is about 
26,000 SHP, whereas at 20 knots the surface 
propulsion requirement has increased to about 
43,000 SHP as against about 29,000 SHP for 


the submarine. 
With regard to ship size it might be expected 


that very large ships would give a large im- 
provement in the competitive position of nu- 
clear propulsion, since the specific capital cost 
(dollars per installed shaft horsepower) of nu- 
clear plants tends to decrease rather rapidly 
with plant output. However, for a given surface 
speed the propulsion power requirement in- 
creases rather slowly with ship displacement, 
and it is probably only when larger displace- 
ment can be accompanied by a substantial speed 
increase (by a factor of 2 or more) that an im- 
portant gain can be made. Actually, reference 1 
concludes that nuclear propulsion is particularly 
attractive in trades that require small ships on 
long trade routes (see the following section). 


Current Construction Projects 
and Design Studies 


There are currently two nonmilitary nuclear 
vessels under construction, the N. S. Savannah’ 
in the United States and the icebreaker Lenin‘ 
in the USSR. Of these, the Lenin is obviously 
not a cargo vessel. The Savannah, designed for 
passenger-cargo service, may be compared with 
other vessels of this type with respect to costs. 
The major characteristics of the reactors for 
these two vessels are listed in Table I-1 along 
with reactor characteristics from a number of 
ship propulsion design studies. 

It has been stated® that an objective of the 
Savannah program is to gain general world 
acceptance of nuclear ships in commerce, and 
consequently the Savannah design includes a 
high standard of safety and reliability. That 
the design is conservative with respect to factors 
which affect power density, coolant circulation 
rate, and system weight is evident by compari- 
son with the icebreaker characteristics in 
Table I-1. 

The Savannah is scheduled for delivery® in 
1960. The ship is a single screw combination 
passenger-cargo vessel, 595'/, ft long, with a 
78-ft beam, a 29'/,-ft draft, and a 21,800-ton 
displacement. It will accommodate 60 to 100 
passengers and carry 10,000 tons of cargo. 
The designed normal speed is 20 knots.°® 

The ship is to be powered by a 74 Mwit) 
pressurized water reactor with two main cool- 
ant loops and steam generators which will pro- 
duce 22,000 SHP at the propeller. The reactor 
vessel will consist of a 26-ft cylindrical shell 
98 in. in inside diameter, with a hemispherical 
head at the bottom and a dished head on top. 
The nuclear plant, including the containment 
and shielding, is located amidship. The ship’s 
conventional propulsion turbines will utilize 
saturated steam. Emergency power for ships 
service, reactor cooling, and take-home power 
will be furnished by two 750-kw diesel gen- 
erators and two auxiliary boilers which bring 
take-home power to 3,000 SHP. 

The primary system is contained in a cylin- 
drical steel vessel, 35 ft in diameter by 51 ft 
long and 2'4 in. thick, lying fore and aft along 
the centerline of the ship. The vessel is de- 
signed to resist the maximum pressure (186 psi) 
that could result from the flashing of all the 
primary water to steam. 
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The primary shielding consists of 33 in. of 
light water and 3 in. of lead aroundthe pressure 
vessel. The secondary shielding surrounds the 
containment vessel and consists of 6'/, in. of 
lead and 6'/ in. of polyethylene covering the top 
half of the vessel, a 4-ft-thick vertical con- 
crete wall around the lower sides of the vessel 
and a 5-ft thickness of salt water ballast and 
feed water tanks under the bottom of the vessel. 

The reactor is protected against collision 
damage by a collision mat consisting of averti- 
cal wall of alternate layers of 1-in.-thick steel 
and 3-in.-thick redwood around the concrete 
shield wall. The mat appears to be between 2 
and 3 ft thick. 

The Savannah is being built for a fixed price 
(conventional practice for shipbuilders) includ- 
ing the ship and nuclear power plant. The total 
capital cost excluding the fuel but including 
some research and development is approxi- 
mately $30,500,000 and is composed as follows:* 


Reactor and power plant $14,500,000 
Shielding and containment 1,500,000 
Hull and outfit 14,500,000 


The design studies listed in Table I-1 cover 
a range of proposed applications as well asa 
range of reactor types. The Japanese, for 
example, have considered the use of a 44,000 
SHP, 20,100-ton gross (29,070 tons displace- 
ment), nuclear-powered vessel to carry emi- 
grants from Japan to South America.* Comparing 
the cost of a PWR powered vessel to a conven- 
tional ship they concluded that, over a 20-year 
period of operation, the financial rate of return 
for the nuclear vessel would be a factor of 2 
lower than that for a conventional vessel (nu- 
clear 17 per cent vs. conventional 34 per cent). 

The British studies’ (Table I-1) report results 
of an economic comparison of a conventional 
with a nuclear-powered supertanker. Although 
power plant capital cost was between 1.5 and 
2 times that required for a conventional ship- 
board power plant, the specific fuel cost (dol- 
lars per SHP hour) was less for nuclear- 
powered vessels in excess of about 25,000 SHP. 
The operating costs for nuclear-powered vessels 
showed that, for fuel oil at acost of 8 cents/gal, 
the nuclear ship broke even (based on operating 
costs) at a cargo deadweight of about 73,000 
tons for a nuclear fuel burn-up of 2700 Mwd/ton. 
Increasing the burn-up to about 5400 Mwd/ton 
decreased the break-even cargo deadweight to 
44,000 tons. 


Atomics International’ has studied the use 
of an organic-moderated reactor (OMR) in a 
20,000 SHP, 38,000 deadweight ton tanker. The 
vessel sails from Kuwait to Philadelphia via 
either the Suez Canal or the Cape of Good Hope. 
The report considers two reactor systems and 
the figures given here and in Table I-1 apply 
to the more advanced, or subsequent, reactor — 
“assuming several other similar maritime OMR 
plants have been constructed so that engineering 
costs are diminished.” The complete vessel cost 
for the OMR powered system was about 25 per 
cent greater than that for a conventional vessel. 
Operating costs were a function of route. The 
nuclear tanker traveling via the Suez Canal 
route operated at a cost about 14per cent higher 
than that of the conventional ship; for the 
voyage around the Cape of Good Hope the op- 
erating costs were about the same. 


American-Standard Study 


By far the most extensive analytical study of 
nuclear ship propulsion has been made by 
American-Standard.'! Four different reactor 
types were studied over ranges of size, appli- 
cation, and trade route variables, and economic 
predictions were made as of the year 1957 and 
for 1965 and 1970. The entire analysis covered 
some 30,000 individual cases. Some of the main 
conclusions are quoted in the following: 


1. Nuclear ships can economically compete now 
with conventional ships on long trade routes at high 
speeds. 

2. The competitive position of nuclear ships will 
improve in 1965 and 1970. 

3. Nuclear fuel costs are, in most instances, 
lower than conventional fuel costs. 

4, Nuclear ship capital costs are 10 to 50% higher 
than capital costs for conventional ships (except for 
gas-cooled reactor powered ships). 

5. Capital costs for nuclear ships will decrease 
with respect to capital costs for conventional ships 
in 1965 and 1970. 

6. Ore carriers and tankers have the best po- 
tential for the economic application of nuclear 
power to merchant ship propulsion. 

7. Nuclear power is particularly attractive in 
trades that require small ships on long trade routes. 

Within the framework of current maritime prac- 
tice, ore-carriers have the greatest area of eco- 
nomic potential for nuclear propulsion. Vaiuable 
cargo tonnage in nuclear ships of this type is not 
sacrificed for fuel space. In addition, this type of 
ship is usually large and requires a large amount 
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of power, which nuclear plants can provide at a 
relatively low fuel cost. 

Ore carriers powered with organic-moderated or 
boiling-water reactors in 1957 can move cargo on 
16,000- to 20,000-mile round trip routes at 19 to 22 
knots at the same unit cargo cost as conventionally 
powered ships. By 1965, large nuclear ore carriers 
operating at 21 to 22 knots on long trade routes will 
be able to carry cargo at about nine-tenths the unit 
cargo cost incurred by counterpart conventional 
ships at the same speed. Moreover, these large nu- 
clear ships will be competitive on a cost per ton 
basis at 18 to 19 knots on the long trade routes, and 
at 19 to 20 knots on routes as short as 12,000 miles 
round trip. Further improvement in the competitive 
position of the nuclear ships occurs in 1970. 

After ore carriers, tankers are the next promis- 
ing application for nuclear power in merchant ships. 
Large nuclear tankers in 1957 can carry cargo on 
20,000-mile round trip routes at lower unit cargo 
costs than conventionally powered tankers, at speeds 
greater than 22 knots. In 1965, the area of competi- 
tiveness for such ships will begin at 12,000-mile 
trade routes. By 1970, nuclear tankers will be more 
economic than conventional ships on 20,000-mile 
trade routes at any speed over 17 knots, and the 
larger nuclear tankers will be economically at- 
tractive on routes as short as 8,000 miles round 
trip. 


Third and fourth in order of economic at- 
tractiveness are nuclear general-cargo ships 
and nuclear passenger ships. 


The particular nuclear propulsion systems which 
offer the earliest opportunity for achieving a com- 
petitive position with fossil-fueled powered plants 
are the organic-moderated and boiling-water re- 
actor systems. In all the cases mentioned above, 
these two types of reactor propulsion systems are 
the ones that are able to compete with conventional 
vessels. Both reactor types present a promising 
potential for the early economic introduction of nu- 
clear technology into the field of merchant ship 
propulsion. From the standpoint of their compara- 
tive developmental status, the most attractive sys- 
tem for nuclear ships in the immediate future is the 
boiling-water reactor, and the most attractive sys- 
tem for the near future is the organic-moderated 
reactor. 


The conclusions summarized above are, onthe 
whole, quite encouraging. There is some un- 
certainty as to the implication of the estimates 
for “present” reactors, for it is stated that the 
reactors are assumed to be fully developed for 
the application considered and that this develop- 
ment would require some further work beyond 
that which has been done to date. Perhaps the 
most informative results for the reactor de- 





signer are those which specify the major per- 
formance and cost characteristics of the reac- 
tors which were arrived at in the study. These 
may be considered as targets which must be 
attained if the economic relationships predicted 
in the study are to hold. Some of these char- 
acteristics are given briefly below. 

In the case of the ore carrier, which was 
judged to be the most promising application, 
the nuclear reactor and machinery cost for a 
boiling water plant of 45,000 SHP was esti- 
mated to be $13,000,000 or $289 per SHP. 
This corresponds to $387 per kw(e) if the 
power output is converted to electrical units. 
This cost applies to a natural circulation 
plant provided with a containment structure 
which would contain all the primary coolant in 
the event of a rupture of the reactor primary 
system. 

Again for the case of the boiling water reac- 
tor in 1957, the fuel cost is estimated to consist 
of a unit variable cost of 0.92 mill/kw(t)-hr 
plus a fixed cost (inventory) of $720 per day 
(for a plant of 45,000 SHP). These estimates 
are based on a fuel lifetime of 10,000 Mwd/ton 
and cylindrical UO, elements with stainless- 
steel jackets. A fuel fabrication cost of $13 per 
pound of uranium (from UO, powder to finished 
elements) is estimated. A plutonium buy-back 
price of $30 per gram is also assumed. The 
plant, operating at 800 psia, gives a thermal 
efficiency of 32.1 per cent. 

At the power level of 45,000 SHP, the operat- 
ing weight of the plant including reactor, con- 
tainment, shielding, and propulsion machinery, 
again for the boiling water reactor, is estimated 
at about 70 pounds per SHP. This figure may be 
compared with that of 365 pounds per SHP for 
the N. S. Savannah. 

With respect to reactor types, the study did 
not show great differences between the pres- 
surized water, boiling water, and organic- 
moderated types; only the gas system, which 
employed a gas turbine, was more expensive 
by a large fraction. All of the reactors studied 
employed partially enriched fuel (except one 
family of pressurized-water reactors which 
was investigated with fully enriched fuel). Both 
the pressurized water and boiling water reac- 
tors used fuel elements composed of steel- 
jacketed UO, pellets. The organic-moderated 
reactor was designed with metal fuel elements, 
jacketed in aluminum. The life expectancy in 
this case was only 3,000 Mwd/ton, but the 
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rather high conversion ratio attained, along 
with the $30 plutonium price, made the fuel 
cycle competitive. The cost of organic make-up 
contributed about 0.12 mill/kw(t)-hr to the 
power cost. The gas-cooled reactor employed 
graphite moderation and steel-jacketed UO, fuel. 
The coolant (and turbine working fluid) was 
helium. It was stated that the induced radio- 
activity in CO, would require shielding of the 
propulsion machinery and would result in a 
weight penalty of 150 to 400 tons over the re- 
actor power level design range of 10 to 200 
thermal megawatts. 


Hazard Considerations 


The hazards peculiar to nuclear systems affect 
the outlook for nuclear ship propulsion in a 
number of ways. The first and most obvious 
effect is the present uncertainty as to the gen- 
eral acceptance of nuclear ships by national 
and local port authorities outside the United 
States. Unless such acceptance is well estab- 
lished the ownership of a nuclear ship will be a 
risky financial venture regardless of other 
economic factors. Presumably the primary use 
of the N. S. Savannah will be to establish prece- 
dents for the acceptance of nuclear-powered 
ships in foreign commerce. 

A second question related to hazards is that 
of the cost and availability of third-party 
liability insurance. In the American-Standard 
study! the total insurance charges were said 
to be 60 to 100 per cent higher for nuclear 
than for conventional ships, and to constitute 
from 12 to 13 per cent of the total daily fixed 
charges. However, it cannot now be said that 
reliable estimates of insurance costs can be 
made. The problem of insurance for nuclear 
vessels is being studied by the Nucleonics 
Committee of the American Institute of Marine 
Underwriters.'° 

Finally, hazards considerations and the degree 
to which precautions are taken affect the capital 
cost of the nuclear ship. The N. S. Savannah, 
with a separate containment vessel and collision 


mats around a relatively low-power-density re- 
actor, probably represents the practical limit 
to which safety precautions can be carried with 
the pressurized-water reactor. These pre- 
cautions undoubtedly add substantially to the 
cost and weight of the plant. It has been stated 
in reference 3 that according to recent design 
studies there is good assurance that ship struc- 
ture can provide adequate protection against the 
release of radioactivity without the use of a 
separate containment vessel, and that this ques- 
tion of the containment feature will become 
pivotal in the development of truly economic 
nuclear ships. 


References 


1. D. L. Conklin et al., Economics of Nuclear and 
Conventional Merchant Ships, ASAE-S-5, June 
1958. (Unclassified AEC report.) 

2.* Michiya Shigemitsu, Nuclear-powered Subma- 
rine Tanker, A/CONF.15/P/1320. 

3.* W. H. Zinn and R. P. Godwin, The Use of Nu- 
clear Energy for Purposes Other Than the Gen- 
eration of Electricity, A/CONF.15/P/1831. 

4.* N. S. Khlopkin, Atomic Icebreaker “Lenin,” A/ 
CONF. 15/P/2140. 

5.* R. P. Godwin and D. L. Worf, Design Considera- 
tions in Nuclear Merchant Ships, A/CONF.15/ 
P/1023. 

6.* S. Takeuchi, T. Okamura, and S. Murakami, 
Nuclear Powered Emigrant Ship, A/CONF.15/ 
P/1319. 

7.* R. V. Moore and C. E. Iliffe, Nuclear Propulsion 
for Ships, A/CONF.15/P/266. 

8.* J. Hainzelin et al., Possibilities and Conditions 
of Use of Water or Gas Type Reactors ina 
Tanker, A/CONF.15/P/1137. 

9. Economic Analysis of an OMR-Powered 38,000- 
DWT Tanker, NAA-SR-1879, March 1957. (Un- 
classified AEC report.) 

10. C. G. Cornwell, Nuclear Powered Ships for 
American Ship Operators, TID-7539, September 
1957. pp. 51-57. (Unclassified AEC report.) 





* Second United Nations International Conference on 
the Peaceful Uses of Atomic Energy, Geneva, Septem- 
ber 1958. (Available from the Office of Technical 
Services, Dept. of Commerce, Washington 25, D. C.) 











REACTOR PHYSICS 





Reactor Physics Constants 


Reference 1 is the first publication of the Reac- 
tor Physics Constants Center, which is located 
at Argonne National Laboratory. The primary 
objective of the Constants Center is stated tobe 
the periodic compilation of the latest and best 
values of the constants, recipes, formulas, etc., 
which are necessary to calculate reactor charac- 
teristics. The compilations, as reported by this 
initial one, which covers data available up to 
Mar. 15, 1958, and some later data, are to be 
quite comprehensive and will undoubtedly be- 
come the first recourse of the reactor design 
physicist. The scope of the work can best be 
summarized by a list of major subdivisions: 
Fission Process, Selected Cross Section Data, 
Constants for Thermal Homogeneous Reactors, 
Lattice Constants for Thermal Heterogeneous 
Systems, Intermediate Reactors, Fast Reactors, 
Shielding Constants, Constants Related to Inter- 
pretation of Experimental Data, and Properties 
of Elements and Reactor Materials.' Obviously 
the magnitude of the task of arriving at “best’’ 
values for the various constants, recipes, and 
formulas requires that this part of the objective 
of the Center must be a long-term one. However, 
the work accomplished to date, as represented 
by this first report, is highly useful, andthe ac- 
tivity of the Center fills an obvious need in the 
reactor program. 


Neutron Slowing Down 


A number of recent papers have given new ex- 
perimental and theoretical information on neu- 
tron slowing down. Most of the information con- 
sists of values of the age for neutrons of various 
initial energy distributions in various modera- 
tors. However, reference 2 is a review of the 
subject of slowing down inhydrogenous modera- 
tors. 


Hydrogenous Moderators 


The pursuit of the classic discrepancy between 
theory and experiment in the slowing down of 


fission neutrons in hydrogenous media is con- 
tinued in reference 2, which summarizes devel- 
opments in this area since the review paper of 
the 1955 Geneva conference.* The following sec- 
tion is a condensation of the first part of that 
report and can give only a rather qualitative pic- 
ture of the subject. Reference 2 and its refer- 
ences (some of which are repeated here) should 
be consulted as a critical review of the subject. 

The unweighted average of three measure- 
ments‘~® of the age to indium resonance energy 
of fission neutrons in light water is 7 = 30.5 + 
1 cm’,* whereas the calculated values are dis- 
tributed as T = 26.0 + 0.5 cm’. Thus the calcula- 
tions consistently result in values that are low by 
about 13 per cent on the average. No reasonable 
permutation of theoretical assumptions and ap- 
proximations seems to have any significant effect 
on this discrepancy, and sources of experimental 
error have been diligently pursuedand carefully 
evaluated. However, some sources of experi- 
mental error have beenonly recently considered 
and no uniform manner of treatment of them has 
yet evolved. These include: 


1. Perturbation by the Source Backing. After 
examining some representative modes of treat- 
ment of this boundary effect, the authors? con- 
clude that the results of the U.S. measurements, 
and probably those reported in the USSR, have 
been adequately corrected in this respect. 


2. Source Perturbations 

Geometric transformations: The authors’ con- 
clude corrections are adequate. 

Displacement of medium by source: The authors’ 
conclude the effect is small. 

Source absorption and scattering: Not signifi- 
cant. 





*An age of 27.65 cm’ was recently measured by 
R. C. Doerner and F. H. Martens (ANL) and reported 
in a paper, Age of Fission Neutrons in Water,at the 
December 1958 meeting of the American Nuclear So- 
ciety, Detroit, Mich. Their results, a summary of 
which is given in Transactions of the American Nu- 
clear Society, 1(2): 30 (December 1958), appeared too 
late for inclusion in this discussion. 
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Source perturbation of background: Only refer- 
ence 4 explicitly states correction. 


3, Detector Perturbations 
Self-shielding of black foils, leading to increased 
proportion of higher resonance activations: 
Corrections appear to be adequate. 
Current-to-flux transformation: Flux age higher 
than measured current age by ~0.5 cm’. 
Finite cadmium cross section at indium reso- 
nance: 


T¢ — Texp ~0.6 cm? 

The authors’ conclude that corrections for 
detector effects are all in the undesired direc- 
tion and barely significant experimentally. 

Sources of error are examined in the several 
“rigorous”* analytical methods that exist for 
computing age from a point isotropic source in 
an infinite homogeneous medium. Truncation and 
round-off errors andstatistical uncertainties (in 
a Monte Carlo calculation) can be made as small 
as desired, and are currently anorder of magni- 
tude smaller than the discrepancy betweentheory 
and experiment. Unknownerrors also arise from 
uncertainties in input data (cross sections and 
spectra of source neutrons). 

The Monte Carlo method has been added to the 
Fourier transform and moments methods pre- 
sented at the 1955 Gevena conference.’ The net 
effect, however, is to leave the 1955 theoretical 
age values unchanged. The transformation for 
slowing down density age to flux age is derived, 
and the detailed assumptions, corrections, and 
input data used by various workers are dis- 
cussed. In particular, values of oxygen total 
cross sections and average cosine of the scat- 
tering angle in the center-of-mass system () 
have been revised’ over the 1953 values® in gen- 
eral use at the time of the 1955 Geneva confer- 
ence. Changes in orbital angular momentum as- 
signments for two of the resonances have radi- 
cally altered the , curve near those resonances. 
These revisions have changed its sign between 
resonances from predominantly negative to pre- 
dominantly slightly positive. The net effect on 
age, however, is small.° 

Further work on oxygen cross sections above 
2.5 Mev is still badly needed; however, acalcu- 





*«‘Rigorous’? here means that no physical assump- 
tions or approximations are made to the Boltzmann 
transport equation. 


lation’® on the basis of aninterim compilation’ of 
oxygen cross sections from 0 to 18 Mev, based 
on all currently available data but with inelastic 
scattering counted as an absorption, has not al- 
tered the 1955 value for the theoretical age. A 
Fourier transform analysis’! using an independ- 
ent evaluation of the oxygen data and including 
inelastic scattering likewise leads to the same 
result. The authors’ conclude that refinements of 
the oxygen data are unlikely to alter substantially 
the calculated age. 


Although the hydrogen scattering cross section 
is definitely anisotropic at 14 Mev, the low- 
energy data are generally interpreted as indi- 
cating isotropy. These data are sufficiently 
rough, however, to be consistent with a 10 per 
cent “anisotropy factor,” 5, where 


_ _9(0) 
~ o(1/2) 





as measured in the center-of-mass system, For 
a symmetrical angular distribution in the center- 
of-mass system, as indicated by the 90-Mev 
data, the average cosine of the angle of scatter- 
ing in the laboratory system, a, is increased by 
0.01276. Furthermore, considerations of nuclear 
forces indicate complete isotropy below 5 Mev. 
The authors’ conclude that the net effect on the 
age of anisotropy of the scattering angular dis- 
tribution must be imperceptible. 

Uncertainties in the fission spectrum are 
shown to be inadequate to account for the dis- 
crepancy between theoretical and experimental 
age. Moments method calculations’ of truncated 
spectra indicate that the discrepancy can be 
eliminated only by neglecting all neutrons below 
600 kev, or about 20 per cent of the total. A 
perturbation approach” shows that the varia- 
tional derivative of the flux age, 7;, with respect 
to the source spectrum, N(E), is 


1 6 _uE)_, 
Ts o NE Tf 


where T(E) is the age from a monoenergetic 
source of energy E; an increase of 60 per cent 
in the spectrum above 2 Mev would be required 
to eliminate the discrepancy. Other types of 
variation lead to even more pessimistic conclu- 
sions. Methods relating the age to the average 
source energy requires 16 per cent increase in 
E, which would be entirely incompatible with 
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available data. The largest defensible increase 
in E would be of the order of 5 per cent, corre- 
sponding to an increase in age of 1.3 cm’. 

In summary, the authors’ conclude that the 
present values for the experimental and theoreti- 
cal age of fission neutrons in water are practi- 
cally equivalent to the 1955 values. The bases 
for the theoretical values are in essence simple 
and clear-cut, and the possible avenues of attack 
seem to have been carefully explored. The most 
uncertain theoretical factor is probably the shape 
of the source spectrum. The most significant 
experimental uncertainty is the absorption of the 
slowing-down neutrons in the source. It seems 
unlikely that the sum of these two effects can 
bridge the gap between theory and experiment. 

The age for neutrons in hydrocarbon modera- 
tors is of interest in connection with the dis- 
crepancy discussed above (as well as intrinsi- 
cally) because the cross-section uncertainties 
for carbon are less than those for oxygen. Carbon 
has no resonances below 2 Mev, and the angular 
distributions are nearly isotropic up to that en- 
ergy. Preliminary measurements” on the age in 
diphenyl indicate an age of 54 +3 cm’, after a 
correction of 8 per cent for transverse leakage. 
The experiments are to be repeated at higher 
power and with a larger tank. A Fourier- 
transform calculation” in the B, approximation 
produced an age of 46.3 cm’, a discrepancy of 
the same per cent magnitude and direction asin 
H,0O. 

A number of age determinations, in addition to 
those listed in reference 1, have been made for 
neutrons from nonfission sources. These data 
also are interesting as additional evidence re- 
lating to the discrepancy between theory andex- 
periment. Reference 2 points out that these data 
have their own peculiar uncertainties, of a mag- 
nitude great enough to weaken greatly the 
strength of any conclusions that might be drawn 
from them. Nevertheless the reference observes 
that the general better agreement between theory 
and experiment for these sources suggests that 
the discrepancy for fission neutrons in water may 
be rooted in some peculiarity of a uranium fis- 
sion source. The data on the nonfission sources 


are summarized below. 

Age measurements have been made at Han- 
ford’ for sodium-beryllium photoneutrons in 
water and kerosene. This source has a spectrum 
of average energy 0.966 Mev, with a full width 
of 3 per cent of the mean energy, and thus ap- 


proximates closely a monoenergetic source. A 
moments calculation’ yielded an age to indium 
resonance of 13.25 cm? in water as comparedto 
a measured value of 13.7 + 0.6 cm’. The age vs. 
source energy curve of Coveyou and Sullivan, 
using the same oxygen data, gives a value of 
14.3 + 0.3 cm’. No calculations have been made 
to check the measured kerosene value of 13.0 + 
0.4 cm’. 


The thermal migration area of the delayed 
neutron spectrum of N"’ was reported by Shure 
and Roys'® as 23.1 cm’. Assuming a migration 
area of 8.7 cm? from indium resonance to thermal 
absorption, the age to indium is 14.4cm? as com- 
pared to theoretical values of 13.6 cm? in the 
Selengut-Goertzel approximation* and 15.2 cm? 
by the “rigorous” method of Coveyou and Sulli- 
van,'® both using the Hayward’® spectrum. Mo- 
ments method calculations are in progress. 


Two experiments were performed at the Na- 
tional Bureau of Standards with neutrons froma 
particle accelerator. One used the 250-kev D-D 
reaction”® to produce 3-Mev neutrons; the other 
used the D-T reaction to produce 14.1-Mev neu- 
trons. The latter approximated a point isotropic 
source, but the oxygen cross sections at that en- 
ergy are so uncertain that the experiment did not 
make possible a fair test of theoretical calcu- 
lations. The measured age was 150+6 cm’, 
Coveyou and Sullivan’ computedan age of 170.4 + 
3.5 cm? at 13.6 Mev counting all nonelastic colli- 
sions in oxygen as absorptions. Using somewhat 
different but unspecified oxygen data, Reier 
et al.'! computed ages of 130 cm? including in- 
elastic, scattering and 185 cm? with zero inelastic 
scattering. 


The D-D experiment”? was complicated by the 
anisotropy of the source in both intensity and 
energy and by the presence of the beam-tube of 
the accelerator. Comparable theoretical results 
do not exist, but a Monte Carlo calculationof the 
flux at 90° to the particle beam, including the 
effect of the evacuated beam tube, is in prog- 
ress.’! This calculation should be directly com- 
parable to the experiment. 


Measurements”*'!»”* of the age to indium reso- 


nance energy of the broad-spectrum Po-a-Be 
source in water give an average age of 56.3 cm’, 
as compared to computed values!! of 53.5 cm? 
and 58.2 + 1.8 cm’, the latter obtained’ by aver- 





*Calculations by R. L. Hellens. 
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aging monoenergetic ages’? over a Murray” 
spectrum. The theoretical values differ by 5per 
cent below and 3 per cent above experiment, re- 
spectively, even though uncertainties in source 
spectrum and oxygen cross sections exist inin- 
tensified form. Measurements of average spec- 
trum energy vary by as muchas 20per cent, yet 
a Fourier transform calculation"! in the P, ap- 
proximation using the MUFT code and four dif- 
ferent spectra led to a variation in age of only 
2 per cent. Furthermore, 50 per cent ofthe Po- 
Be neutrons are above 4 Mev in energy, where 
knowledge of oxygen cross sections is sketchy, 
whereas only 10 per cent of fission neutrons ex- 
ceed this energy. In addition, the manner of 
treating disappearance processes and inelastic 
collisions led to a 10 per cent variation in 7, 
whereas the use of two different sets of oxygen 
data produced’ a 9 to 11 per cent variation in r 
at 3.6 and 4.7 Mev. In the face of uncertainties 
of this magnitude, it cannot be claimed that there 
is any discrepancy between theory and experi- 
ment for Po-Be neutrons in water. 

The measured’ age of 148 cm? for Po-Be neu- 
trons in 99.4 per cent D,O compares well with 
the calculated'® age of 150 cm’, although the as- 
sumption of isotropic scattering in the center- 
of-mass system at all energies suggest chance 
agreement. 

Measured values"! of Po-Be neutrons in water- 
metal mixtures agree well with MUFT values 
except for zirconium-water mixtures, where the 
knowledge of zirconium inelastic cross sections 
is poor. Incidentally, the measurements indicate 
that volume for volume, iron is almost as good 
a moderator as water, due to inelastic scat- 
tering.’ 

The neutron spectrum ofthe Po-a-B source is 
poorly known but is considerably narrower than 
that of Po-a-Be.* The theoretical age in water”4 
for various spectra ranges from 30 to 36 cm’, 
vis-a-vis a measured’ age“ of 39.7 + 1.2 cm’, 
while that in a 2:1 Al-H,O lattice”® varies from 
109 cm? to 123 cm? vs. a measured age of 130 + 
3.9 cm’. Thus a discrepancy of the same order of 
magnitude and sense exists for the Po-B-H,O 
situation as for the fission case. 


D2O-H20 Mixtures 


The age in mixtures of D,O andH,Ois of con- 
siderable practical importance in the design of 
D,O-moderated reactors where the possibility of 
contamination of the D,O by H,O exists. Meas- 


urements have been reported for mixtures con- 
taining from 99.8 per cent D,O” to as low as 
48.6 per cent D,O.* For the mixtures containing 
small percentages of H,O the age is very sensi- 
tive to the H,O content, decreasing by about 4 
cm’ for 1 per cent increase inH,O content. When 
the H,O content becomes large the age changes 
only slowly with composition. By the time the 
H,O content has been increased to 51.4 per cent 
the age to indium resonance has dropped froma 
value of about 109 cm? (for 99.8 per cent D,O) to 
38.6 cm’, which is higher by only 8 cm? than the 
value for pure H,O. 

Existing theoretical calculations’~"* of the age 
for pure D,O give results somewhat higher than 
measured values, probably because the prefer- 
entially backward scattering by deuterium was 
not taken into account. The calculated values for 
D,O-H,O mixtures are inagreement with experi- 
ment at about 60 per cent D,O content, and, of 
course, fall below experiment for pure H,O. 
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Graphite and Graphite-Bismuth Mixtures 


The age to indium resonance for fission neu- 
trons has been reported, 7° from exponential ex- 
periments, as 285.33 + 0.4 cm? for graphite of 
density 1.674 g/cm*. This would correspond to 
a value of 311 cm’ for graphite of density 1.60 
g/cm® and may be compared to the value from 
the Reactor Handbook”’ of 312.5 +0.5 cm? for 
graphite of that density. A value of 449.8 + 2.0 
cm’ has been reported'’ for a graphite-bismuth 
stack, having a graphite-bismuth volume ratio of 
2. The graphite density was 1.674g/cm’, andthe 
stack contained some aluminum in the form ofa 
foil wrapping for each bismuth bar. This meas- 
urement is useful in connection with the design 
of liquid-metal fuel reactors. 


Monoenergetic Neutrons in lron 


and Stainless Steel 


Roos”® has recalculated Horton’s*® ages to 
thermal with recent cross-section data®’ and 
finds that the sharp decrease in cross section at 
27 kev in pureiron leads to increases by a factor 
of 2 in the age to thermal for source energies in 
excess of 27 kev. The factor could well become 
larger as experimental resolution of the cross- 
section “window” improves. In stainless steel, 
the window is partly covered by the alloying ele- 
ments, and the changes over Horton’s values are 
much less drastic. 
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Physics of Low-enrichment 
Light-water Reactors 


The physics of the light-water-moderated and 
-cooled power reactor has been the subject of 
rather extensive experiment and analysis. Al- 
though-a considerable amount of research and 
development has been on a prototype basis and, 
as such, the information does not have general 
application, there is at present a large body of 
experimental data available against which meth- 
ods of analysis can be tried. 


Recent papers have summarized the status of 
theory and experiment for light-water lattices*® 
and the application of digital-computer tech- 
niques to the calculation of lighi-water-moder- 
ated reactors.*' A digest of this material is be- 
yond the scope of the present review, but the 
following section will serve as a guide to the 
current sources of experimentai and theoretical 
information. 


Experimental 


During the past few years, experimental tech- 
niques have been greatly improved, and rather 
difficult measurements are now made inroutine 
fashion. Further, there has been a greater em- 
phasis on the performance of experiments that 
lend themselves directly to prescribed analytical 
methods. Much of the experimentation on H,O- 
moderated lattices has been done at Brookhaven 
National Laboratory (BNL) and at the Westing- 
house Atomic Power Division (Bettis). 


The BNL experimental measurements on 
metallic rod lattices have been described in 
references 31, 32, and 33. Fuel rods with0.250, 
0.387, and 0.600 in. diameters, clad with 9.028- 
in. 2S aluminum, with a nominal 0.005-in. air 
gap between fuel and clad, have been positioned 
such that the water-uranium volumes ranged 
from 1 to 4. The U?** enrichments studies (by 
weight) were 1.027, 1.143, and 1.299 per cent. 


The experimental data at Bettis have beenob- 
tained by critical assembly techniques,*4~*® 
Measurements were performed on the identical 
metallic rods used in the BNL experiments for 
the 0.387- and 0.600-in. cases with 1.143 and 
1.299 per cent enrichments. Also, oxide meas- 
urements were made for 1.311 per cent enrich- 
ment with UO, density 7.52 g/cm® and 10.53 
g/cm*®. The oxide fuel diameter was 0.385 and 


0.601 in., and the water to oxide volume ratio 
ranged from 1 to 2.5. 

For all the lattices described above, buckling 
measurements have been made, as wellas deter- 
minations of lattice parameters, including fast 
effect, resonance escape probability and disad- 
vantage factors. In addition, measurements of 
reflector savings, migration areas, and tem- 
perature coefficients were made inthe tempera- 
ture range 20 to 50°C. 

A set of measurements was also made at BNL*® 
on lattices whose fuel was 1 per cent plutonium, 
0.3 per cent U**5, and 98.7 per cent U***. The 
rods were 0.387 in. in diameter and were clad 
like the slightly enriched lattices mentioned 
above. 

Zero-power experiments have been performed 
on the EBWR reactor.*’ These experiments in- 
clude four different plate type uranium-metal 
fuel elements of two-plate thicknesses and in two 
enrichments. One plate is 0.214in. thick, andthe 
other is 0.279 in. thick; bothare clad with 0.020- 
in. Zircaloy-2. The two enrichments are natural 
uranium and 1.44per cent enricheduranium. The 
measurements reported include critical mass 
determinations, control rod calibrations, void 
and temperature coefficient measurements, de- 
terminations of reactivity compensated by steam 
voids, and the effect of boric acid and control 
rods on power generation. 

Extensive measurements have been made on 
the Shippingport reactor’ utilizing zero-power 
mock-up facilities.*°-*° The geometrical con- 
figuration consisted of a seed and blanket. The 
seed was composed of zirconium plates as- 
sembled into subassemblies approximating the 
mechanical design of the PWR seed. The seed 
fuel, in the form of enriched uranium oxide, was 
dispersed in plastic tape and cemented to the 
zirconium plates. The blanket material consisted 
of natural uranium oxide pellets, identical to 
those used in PWR, inserted into aluminum tubes 
having approximately the dimensions of PWR 
blanket fuel assemblies. The principal measure- 
ments in the mock-up facility comprised deter- 
minations of critical configurations, excess and 
shutdown reactivities, power distributions, and 
temperature coefficients. 

Zero-power experiments on a mock-up of the 
Yankee Atomic Power Plant‘! have recently been 
performed.**-‘4 The fuel rods consisted of 
0.300-in.-diameter UO, pellets sealed in type 
304 stainless-steel tubing with a wall thickness 
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of 0.016 in. The U?®* enrichment is 2.7 wt.%. 
The fuel rod spacings were varied to give water- 
to-fuel volume ratios from 2.23 to 3.81. Meas- 
urements of critical size, buckling, reflector 
savings, migration area, disadvantage factor, 
resonance escape probability, and temperature 
coefficient were made for each core. 

Experiments similar to the Yankee plant have 
been performed onthe Savannah Merchant Ship.” 
The fuel consists of UO,, enriched to 4 per cent 
U5 in the form of 0.500-in.-diameter rods with 
a 0.028-in. stainless-steel clad. The fuel has 
been arranged in critical arrays with water to 
UO, volume ratios® of 1.15 and 1.54. 

As a part of the Russian program on pres- 
surized-water reactors*> measurements of lat- 
tice parameters and bucklings for metallic rods 
of approximately 1.35 in. in diameter have been 
made for water to uranium volume ratios*® 
ranging between 0.4 and 2.0. For the latter ex- 
periments spiked and seed-blanket arrange- 
ments were also measured. Measurements of 
the neutron spectrum in uranium-water lattices 
were also made for the 1.35-in. rods.*’ Appar- 
ently these experiments, with large metallic 
rods, were made for general information; the 
Russian pressurized water designs employ oxide 
elements of relatively small diameter. 

Experiments have been performed in the 
United States on several enriched uranium- 
thorium-water critical assemblies*® in connec- 
tion with the Consolidated Edison Reactor design. 
These experiments have included both plate type 
elements containing thorium metal and pin type 
elements containing thorium oxide. Lattice pa- 
rameters, critical masses, and bucklings and 
migration areas were determined for bothtypes 
of assemblies for several uranium-to-thorium 
ratios. 


Organic materials such as diphenyl have mod- 
erating properties quite similar to those of light 
water reduced to an equal hydrogen density. Ex- 
ponential experiments have been performed with 
lattices consisting of 1l-in.-diameter metallic 
uranium rods in liquid diphenyl.*? The fuel used 
consisted of both natural uranium and0,912 at.% 
u**>, For each enrichment, lattices of moderator- 
to-fuel volume ratios ranging from 1.5to4were 
studied. Measurement of bucklings, age, and 
diffusion lengths were made. 


In addition to the lattice experiments de- 
scribed above, the recent revised cross-section 
compilation 


*0 is of course pertinent to the sub- 


ject. Further, measurements in the fields of 
neutron spectra, reactor control, and reactor 
kinetics are also important. Much of the recent 
information on these subjects is summarized in 
the Argonne National Laboratory Reactor Phys- 
ics Constants compilation. ' 


Theoretical 


The “classical” formulation of the theory for 
heterogeneous thermal reactors was developed 
for well thermalized lattices. Slightly enriched 
H,O-moderated lattices, of compositions which 
can attain criticality, are in general not well 
thermalized; i.e., an important fraction of the 
fissions is caused by neutrons of epithermal en- 
ergy, and the thermal-neutron energy spectrum 
is determined by the ratio of absorption to mod- 
erating power as well as by the temperature of 
the moderator. Other peculiarities of H,O- 
moderated lattices are the relatively important 
fast-fission effect, involving interactions be- 
tween fuel lumps, and the non-Gaussian spatial 
distribution of the slowing-down density. Asa 
result of some of the peculiarities, the theoreti- 
cal bases for expecting the approximations in- 
volved in the “four-factor” formulation to be 
good ones no longer exist (although reference 
33 states “It may however be stated that a previ- 
ous conclusion [see reference 32] still seems 
valid: A four-factor interpretation is self- 
consistent if the values of 7 which are used are 
those based on a choice of 1.309 for natural ura- 
nium”). Consequently the current theories for 
H,O-moderated lattices do not involve the four- 
factor approximation. Reference 33 gives a re- 
view of the basic approach used in most of the 
current calculational schemes. Salient points of 
the approach are 

1. The theoretical treatment begins from a 
spherical P, approximation to the transport 
equation, which deals with homogenized neutron 
properties and contains self-shielding factors to 
allow for those features which homogenization 
would ordinarily suppress (higher approxima- 
tions to the transport equation may of course be 
used to calculate the self-shielding factors). 

2. The energy distribution of the neutrons is 
treated by few-group theory——a multigroup 
theory comprising a relatively small number of 
groups. Fission and absorption are allowed in 
all groups; hence, the limitations of the four- 
factor formulation, which conceptually restrict 
the processes of fast fission, resonance absorp- 





REACTOR PHYSICS 13 


tion, and thermal absorption to specific energy 
regions, is avoided. It is stated that often two, 
three, or four energy groups appear to be ade- 
quate for most design purposes. For suchbroad 
groups cross coupling of the groups can be neg- 
lected. The few-group constants may be arrived 
at by solving a simplified problem by many- 
group theory and choosing the few-group con- 
stants to be consistent with the resulting energy 
spectrum. Other papers on the computation of 
group constants are references 52, 53, and 54. 

3. Thermal neutrons are treated by one-group 
diffusion theory; the top energy of the thermal 
group is chosen sufficiently high that no scatter- 
ing out of the group occurs. The thermal-neutron 
energy spectrum is determined bya calculation, 
such as the Wigner-Wilkins, which takes into 
account the equilibrium between absorption and 
moderation. Usually this calculation is made by 
machine, using a code suchas SOFOCATE, which 
has the advantages over the Wigner-Wilkins 
method of allowing absorption cross section 
variations other than 1/v, and of accounting for 
neutron leakage. 


Reference 85 gives a rather detailed descrip- 
tion of the set-up of afive-group model for H,O- 
moderated lattices, the five-group model being 
derived from a 22-group analysis. The modelis 
applied to several lattices for which experi- 
mental data are available, and comparisons are 
made with experimental results. 

Although the basic theory of the reactor is of 
fundamental importance, much ofthe theoretical 
and computational effort on H,O reactors has 
gone into the solution of practical design prob- 
lems which owe their difficulty to their com- 
plexity. Typical problems are the prediction of 
the variation of the spatial power distribution 
over the life of the core, the predictionof reac- 
tivity lifetime when the fuel burn-up is nonuni- 
form, and the calculation of temperature and 
void coefficients of reactivity for nonuniform 
cores. Many of these problems have been at- 
tacked by machine computation. A review of the 
role of the digital computer in the design of 
water reactors is given in reference 51. The 
codes for predicting reactivity behavior solve the 
multi-group diffusion equations. Many one- 
dimensional diffusion codes are now in use. 
These codes include WANDA,» PROD,” and 
GNU.” CANDLE® and BOX are one-dimensional 


codes that include composition changes due to 
the burn-up and buildup of different isotopes. 


Two two-dimensional codes that are now in op- 
eration are CURE™ and PDQ.” These codes also 
include provisions to calculate composition 
changes as a function of lifetime. TRIXY® is a 
three-dimensional code now in operation. TKO 
and FLAME are three-dimensional codes which 
are now being developed. 


Doppler Coefficient for U8 


In reactors which utilize low-enrichment ura- 
nium as fuel, the temperature coefficient of the 
effective resonance integral for U**® is a pre- 
dominant factor in the prompt temperature coef- 
ficient of reactivity. It is important in safety 
evaluations, and it contributes to the reactivity 
change between the low-power low-temperature 
condition and the hot operating condition. Meas- 
urements of the variation in the U’*® effective 
resonance integral with fuel temperature have 
been carried out by several investigators. The 
results of the different measurements are not in 
good agreement relative to the accuracy required 
in reactivity calculations. High accuracy is es- 
pecially important in calculations for reactors 
utilizing natural or near-natural UO, as the fuel 
because ofthe marginal excess reactivity availa- 
ble in such reactors and the high operating tem- 
perature of the UO,.. However, the order of mag- 
nitude of the effect is well established by the 
various measurements and also by caiculations 
based on experimental data on the resonance 
structure of U2°*. The values of dA/AdT inter- 
preted from the various measurements range 
from ~1.3 x 107! per degree centigrade to ~2.1 x 
107 per degree centigrade (between room tem- 
perature and about 650°C), where A is the 
volume-absorption term in the formulation 


A ee) 


RI == + 
: F B mass 


for the resonance integral. In this formulation, 
A and B are constants of the material, and F is 
the resonance disadvantage factor. 

The results of a comprehensive set of meas- 
urements have recently been reported on the 
temperature coefficient of the effective reso- 
nance integral for u’*® over the temperature 
range 20 to 600°C for the metal and 20 to 1000°C 
for UO,. The surface-to-mass dependence of the 
temperature coefficient was investigated by the 
use of fuel slugs of three different diameters, 
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ranging from 8 to 28 mm for both the metal and 
the oxide fuel. Thus it was possible, from analy- 
sis of the data, to determine the separate tem- 
perature coefficients of the two absorptionterms 
comprising the foregoing formulation for the 
resonance integral. The rods were irradiatedin 
a cadmium-covered oven in the reflector of the 
Swedish D,O reactor R1. Measurements were 
carried out by an activation technique based on 
obtaining the intensity, for solutions of the ir- 
radiated samples in nitric acid, of the 105-kev 
X-ray line in Pu*® following the Np?*’ decay. 

In the interpretation of the activation meas- 
urements, the temperature dependence of the ef- 
fective resonance integral was taken to be of the 
form 


‘x | > 


RI(t) = RI(1/v) + I! +a(t- 20)| 


+[1 +8 (t- 20] BS 


The percentage increase in specific activity of 
the irradiated samples was assumed to be line- 
arly dependent on the temperature over the tem- 
perature range investigated. Within the data 
scatter and the estimated probable errors for 
each data point, this assumption is indicated by 
the plotted data to be about as reasonable as any. 
The values given by Hellstrand® for A (7.8 
barns for metal and 11.6 b for oxide) andB 
(30.2 b g/cm’ for metal and 22.8 b g/cm? for 
oxide) and the value for RI(1/v)(1.lb) were used. 
On these bases, the values of @ and 8 were ob- 
tained by a least squares fit tothe reduced data. 
The values obtained for metal are 


a = 0.6 x 10~! per °C 
and 

B = 2.6 x 1074 per °C 
The values obtained for oxide are 

a = 0.1 x 1074 per °C 


B = 3.1 x 1074 per °C 


As interpreted, the experimental results show 
a temperature dependence whichis mainly due to 
the effect of the surface absorption term. 

The temperature coefficients of the resonance 
integral obtained in these measurements are 


compared for U?** metal, on a common basis, 
with the results of other investigators. The re- 
ported results exceed that given by Rodeback®™ 
but are lower than the results of Pearce and 
Walker™ and of Davis.™ 


Effect of Holes on Slowing Down 
and Diffusion Areas 


The presence of voids withina reactor medium 
increases the neutron age and thermal diffusion 
area in the medium and hence results in in- 
creased neutron leakage from the reactor. Fora 
uniform distribution of voids whose dimensions 
are small compared to the neutron mean free 
path in the solid medium surrounding the voids, 
the increased age andthermal-diffusion area can 
be simply obtained from the age and thermal- 
diffusion area values when no voids are present. 
The expression is 


L? = L? (1 + 9)’ 


with a similar expression for r. Here 7 and L” 
are the age and thermal-diffusion area, respec- 
tively, with subscript o indicating the case of no 
voids. The ratio of void volume to solid-material 
volume is 9. 


As the dimensions of the individual void vol- 
umes are increased and become large compared 
with the mean free pathinthe solid material, the 
foregoing simple treatment is no longer correct. 
In this case, the increase in 7 and L? due to the 
introduction of the voids can be substantially 
larger than that indicated by the simple expres- 
sion. Basically, the reason for the greater in- 
crease is that the neutron transport properties 
are governed by the root-mean-square of the 
neutron flights rather than by the arithmetic 
mean free path in the medium. As a further re- 
sult the size, shape, and orientation of the indi- 
vidual void volumes are important considera- 
tions. 


A theoretical analysis of the increase in the 
migration area of neutrons in a reactor by the 
presence of voids was made by Behrens” and 
includes the effect of size and shape of the voids. 
For void geometries characteristic of the holes 
passing through the reactor medium for the flow 
of the coolant in a gas-cooled power reactor, 
(the most important practical case), Behrens’ 
results are expressible as 
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L2 

TE n1+2 0+ Se 
(9) t 


where 


L?, and L’ =the thermal diffusion areas paral- 
lel and perpendicular, respec- 
tively, to the gas-flow passages 

A, = the thermal neutron mean free path 
in the solid medium 

r=the hydraulic radius of the hole 
(twice the volume of the void di- 
vided by the surface area of the 
void) 

Q =a function of the shape of the hole 


A similar expression, with A; equal to the neu- 
tron mean free path during the slowing-down 
process, is appropriate for tT. The factor Q has 
been evaluated by Behrens” for various hole 
geometries. 

Experimental results have recently been re- 
ported on the anisotropic diffusion lengths for 
thermal neutrons in large graphite assemblies 
containing void channels which simulate the gas- 
flow passages in large gas-cooled power reac- 
tors. Prior to operation of the G1 reactor (natural 
U-metal, graphite- moderated)” the reactor was 
partially loaded with fuel, the load being confined 
to one end of the graphite stacking and being 
sufficiently extensive to achieve criticality. The 
thermal-neutron diffusion was observed in the 
nonloaded end. The graphite stacking was pierced 
by gas-flow channels of 7 cm in diameter forming 
a square lattice of cells of 20-cm pitch. Meas- 
urements were also made with the channels filled 
with graphite of the same quality used in the 
stacking. The experimental results and a com- 
parison with Behrens’ theory are 


@ = 0.106 
(1 + o)? = 1.223 


L, = 52.7+ 0.4 cm 


Behrens’ 
Experiment theory 
(Lig, /L* 1.504 + 0.040 1.494 
(L,/L,)” 1.393 + 0.045 1.353 
e=(L,,/L,) 1.08 + 0.03 1.104 





An experimental test of Behrens’ theory is also 
provided by reference 71 in which a large value 
of @ = 0.333 was used. A graphite stacking ap- 
oroximately 66 by 65 by 73 in. was used, and the 
source of thermal neutrons was obtained froma 
graphite pedestal with an antimony-beryllium 
source. Measurements of thermal neutron flux 
were made with and without the void channels (of 
square cross section) running the length of the 
pile. The results obtained and comparison with 
Behrens’ theory are tabulated below. 


o = 0.333 


L, = 54.7 + 1.8 cm 


Experiment Behrens’ theory 
L,, (cm) 94.0 + 3.0 94.6 + 3.1 
L, (cm) 82.2 + 2.5 83.5 + 2.7 
€ 1.305 + 0.03 1.284 


The two separate sets of measurements pro- 
vide, for the first time, experimental confirma- 
tion of Behrens’ analytical results, at least for 
assemblies typical of the large gas-cooled 
graphite- or D,O-moderated reactors fueled with 
natural or near-natural uranium. These results 
indicate that for the relatively large assemblies 
used, the neutron-streaming “end effects” are 
small. Near the ends of the reactor assembly, a 
significant fraction of the neutron flights into the 
gas passages lie in a range of solid angles such 
that the neutron is not again intercepted by the 
reactor medium; since Behrens’ analysis is for 
an infinite medium this end effect is not accounted 
for in his theory. 


In reference 72 the statement is made that 
“the migration area asymmetry associated with 
the presence of large gas-cooling channels (for 
the Calder Hall reactor) was consistent with the 
Windscale observations and greater than the 
simple predictions of Behrens.’’ It was later 
stated that Behrens’ values were used, and 0.2 
per cent reactivity was allowed for the under- 
estimate of neutron leakage from the Calder Hall 
reactor. Since the leakage from Calder Hall is 
only of the order of 3 to 4 per cent, the 0.2 per 
cent reactivity allowance will take care of an 
underestimation of neutron migration area of 5 
to 7 per cent. This is a somewhat larger error 
than is indicated by the results of references 70 
and 71. 
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Recent Measurements 
on DIMPLE and GLEEP 


The experimental program being carried out 
at Harwell using the pile oscillator facilities of 
DIMPLE and GLEEP is described in reference 
73. The work includes measurements of (1) 
thermal cross sections of various purely absorb- 
ing materials of interest, including some of the 
important stable fission-product poisons; (2) 
resonance absorption integrals, of which that for 
Pu“° is of special interest; (3) 1 values for U?* 
and Pu?’ relative to that for U*; and (4) reac- 
tivity changes of irradiated fuel from the 
Windscale and NRX reactors. The results ob- 
tained to date are reported and compared with 
existing information.’? Because of the impor- 
tance of such basic data in reactor calculations 
and particularly of the measurements made in 
several reactor spectra, a summary of the re- 
sults is given here in some detail. 


In natural or near-natural uranium systems it 
is convenient and appropriate to use the cross- 
section convention and pile-spectrum character- 
ization introduced by Westcott in reference 74. A 
1958 revision of reference 74 has recently been 
published by Westcott."* In this formulation, the 
effective neutron cross section is so defined that 
its product with the “2200 meters per second 
flux’’ (i.e., flux given by the total neutron den- 
sity, including the epithermal neutrons, multi- 
plied by v, = 2200 meters per second) gives the 
correct reaction rate directly. Obviously, such 
a formulation is useful only if the epithermal 
contribution is not large; this is usually true in 
natural uranium or near-natural uranium sys- 
tems. On this basis, the eff2ctive cross section 
6 is given by 


6 = 0, (g +rs) 


where 6, =2200 meters per second cross- 
section value 
g =a measure of the departure of the 
cross-section variation from the 1/v 
behavior in the thermal region (fora 
1/v absorber, g = 1) 


4T Or 


st 
Vailg Uo 


with T the neutrontemperature, T, = 
293°K, and o, the resonance integral 





excess over 1/v (for a1/vabsorber, 
s = 0) 

r =a measure of the ratio of the epi- 
thermal to thermal neutron densities 


In this formulation, epithermal flux is assumed 
to be distributed in energy as 1/E; hence, again, 
the applicability only to systems of low epi- 
thermal component is apparent. 

From measurements in a pure Maxwellian 
spectrum (r ~ 0), g 0, is obtained directly. From 
further measurements in a neutron spectrum 
with a known epithermal component (as given by 
the value of r), o, can be deduced from the dif- 
ference between the two results. 

In DIMPLE, the spectrum at the sample posi- 
tion is extremely well thermalized (r < 0.001). 
In GLEEP, r = 0.045 at the sample position. 

Values of o, obtained from the measurements 
for various purely absorbing materials compare 
well with BNL-325 (and its Supplement No. 1) 
tabulation. For the important stable fission- 
product Sm"*, however, the oscillator value is 
42,100 b compared to 66,000 b as givenin BNL- 
325. The 42,100 b value is said to be in good 
agreement with that recently measured™ by a 
different technique. For N,”*", g o, was deter- 
mined as 169 b. Since Np’*’ has been shown to 
follow a near 1/v law, the measurement gives a 
0, value in good agreement with other investiga- 
tors. A value of o, for Pu*4° was deduced as 
370 + 40 b from measurements on a sample 
whose isotropic composition of Pu**, Pu?4?, 
Pu"! and Pu’? was determined by mass spec- 
trometer analysis; comparison is made in the 
reference with the results of various investiga- 
tors. 

Some of the data on resonance absorption inte- 
grals obtained in a hardened neutron spectrum 
in DIMPLE (r = 0.098 obtained by loading of 1.2 
per cent U?*> fuel around oscillator tube) are in 
conflict with those of Macklin and Pomerance.” 
For Hf, o, = 2800 b was measured as compared 
to the 1250 b value reportedin reference 77. For 
Pu“? a value of o, = 11,300 + 1000 b was meas- 
ured, and this value is compared in the refer- 
ence with that of other investigators. The USSR 
parameters list o, (Pu’**) = 10,000 + 1500 b. 

Measurements of relative values of 7 of the 
fissile nuclei have also been made by the pile 
oscillator technique. The oscillator measure- 
ment (by comparison with a standard sample) 
gives Ors, which is related to the absorption 
cross-section o, and the neutron yield 7 by 
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O off = Tq (1 — Wn) 


where W is the ratio of the fission-neutron to 
thermal-neutron importance at the sample posi- 
tion. Because of the similarity of the fission- 
neutron spectra of the fissionable isotopes, Wis 
taken to be the same for U’*>, U?%3, and Pu’9?, 
Hence 


nx (1-52) 
Oa 


where Oj ¢- is given by the pile-oscillator meas- 
urements and where appropriate values of o, 
were obtained from Westcott’s compilation. “4 It 
is stated in the reference that uncertainties in 
0, introduce an uncertainty of about 1.5 per cent 
in the determination of the relative 7 values. The 
results are given in the following table with 
Westcott’s predicted ratios for these spectrum 
conditions shown in parentheses. 

n(Pu**)/n(U*) 


Reactor spectrum n(U*33) /n(U*5) 


DIMPLE 
(~Maxwellian) 1.057 + 0.015 0.974 + 0.015 
(1.103) (0.987) 
GLEEP 1.091 + 0.015 0.980 + 0.015 
(1.102) (0.962) 


For U’*?, the situation is confusing. Onthe basis 
that n (u?*5) in the DIMPLE spectrum is expected 
to be not more than 0.5 per cent higher than in 
GLEEP, while 7(U**) should have the same value 
in both reactor spectra, the DIMPLE measure- 
ments for the 7 ratio should be only about 0.5 
per cent lower than the GLEEP measurement. 
The GLEEP measurement gives reasonably good 
agreement with Westcott’s predicted ratio; how- 
ever, it is stated in the paper® that the use of 
“better data” for o, (U’**) in the analysis of the 
oscillator measurements would lead to a lower 
n ratio. The DIMPLE measurement gives an 7 
ratio some 4 per cent lower than Westcott’s 
predicted ratio. 

The over-all absolute discrepancies in the 
value of n(U’**) are large. If an average of the 
DIMPLE and GLEEP measured values of n(U”**) / 
n(U?**) is taken and if 7(U?*5) is taken as 2.07, 
then 7(U?**) = 2.22. This compares to a value of 
2.28 obtained when using Westcott’s value for 
the ratio. However, the situationis further com- 
plicated by the suggestion in reference 78 that 
n(U?*5) = 2.04, made on the basis of current in- 


terpretations of exponential experiments at 
Harwell; in this case, when the average of the 
DIMPLE and the GLEEP measured ratios is 
used, 7,(U’**) becomes equal to 2.19. Hence, the 
present uncertainty in n(U") is of the order of 
4 per cent. 

The situation with respect to Pu*’’ is con- 
siderably better. The DIMPLE value of n(Pu?*’)/ 
n(U?*5) is in fair agreement with Westcott’s 
recommended 0.987 and agrees well with the 
value of 0.97 given by Harvey and Sanders.” 

Reference 1 (Sec. 2) presents a detailed com- 
parison of 7 values determined by various in- 
vestigations for uv U3, and Pu’: in the 
comparison, plots are given of n(E) vs. neutron 
energy E. 

In a panel discussion at the winter meeting of 
the American Nuclear Society, Weinberg” called 
attention to the great significance of then of U?*? 
for the future of U’** breeding, and stated that a 
concerted effort is being undertaken inthe United 
States to reduce the uncertainty inthis constant. 

GLEEP measurements have also been made 
for a variety of irradiated fuel materials. The 
results are being interpreted using a fairly 
simple model of neutron behavior for the pur- 
pose of deriving reasonable combinations of the 
pertinent nuclear constants with sufficient ac- 
curacy for predictions of long-term reactivity 
changes. A typical analysis is given inthe paper ”® 
for a set of some 20 NRX uranium bars irradi- 
ated to a maximum of ~2400 Mwd/ton. In the 
analysis of the pile oscillator experiments, 
which measure the change in effective cross 
section of the bar with irradiation exposure, it 
is necessary to account accurately for: 

1. The change in thermal neutron blackness of 
the bar with irradiation exposure. 

2. The relative blackness of the unirradiated 
bar and the standard boron sample used. 

3. The growth of Sm“* from the decay of Pm'*® 
after reactor shutdown. 

4. The over-all contribution from the fission 
products. 


239 


Calculations were also made of the changes in 
the isotopic composition of the fuel with irradia- 
tion exposure using an NRX conversion ratio of 
0.76 and using effective cross sections appro- 
priate to the spectrum conditions inGLEEP. The 
heavy isotope cross-section data used in the cal- 
culations were those presented in references 81 
and 82. In order to convert the calculational re- 
sults to reactivity contribution as interpreted 
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from the GLEEP measurements, a calculated 
value of W (ratio of fission to thermal neutron 
importance at the sample position) was used. 

The predicted reactivity variation is in fair 
agreement with the measured variation. It 
should be pointed out that the 7(U’*>) for 2200 
meters per second neutrons was taken as 2.04 
in the calculations (as given by the source of 
cross-section data used). If a value of n(U***) of 
2.07 rather than 2.04 were used, the predicted 
curve shows appreciably less reactivity than 
the measured curve. This illustrates the sensi- 
tivity of the prediction to small variations inthe 
nuclear data. 


Estimated Reactivity Lifetime 
in D,O Reactors 


From a critical study of the changes in iso- 
topic composition and nuclear parameters of the 
fuel with irradiation exposure, it is concluded in 
reference 83 that there is a good possibility that 
fuel exposures greater than 8000 Mwd/tonne* can 
be achieved in a large D,O-moderated natural- 
uranium power reactor using fuel in the oxide 
form. This conclusion is pertinent to any evalua- 
tion of the D,O-moderated reactor as a typeand 
the methods used in arriving at the estimateare 
of general interest in relation to the prediction 
of reactivity lifetime. Obviously, the conclusion 
must be based on an estimate of the neutron 
economy attainable in a practical power reactor. 
The exposure prediction is based also onthe use 
of the bi-directional fuel-push-through refueling 
scheme. The neutron balance, as it varies with 
fuel irradiation, is extremely sensitive to (1) the 
relative rates of depletion of U5 and ingrowth 
of Pu2®®, Pu24?, pu24! and fission products (which 
depend on the effective cross sections of the 
isotopes and the rate of conversion of U*** to 
Pu?%*), and (2) the 7 values for the fissionable 
isotopes. The effective cross sections and 7 
values depend on the neutron spectrum seen by 
the “fuel” during the irradiation exposure. Un- 
certainties in the basic nuclear data andin prop- 
erly accounting for the neutron-spectrum effects 
may introduce rather large uncertainties in pre- 
dicted fuel exposures as limited by reactivity 


considerations. 
In reference 83, the evidence is summarized 


on the basis of which the prediction of a yield 


*The ‘‘tonne’’ is the metric ton. 





greater than 8000 Mwd/tonne is made. One im- 
portant piece of evidence is the agreement ob- 
tained between the results of calculations by 
Kushneriuk and the results of analysis of NRX 
rod No. 683. This natural uranium metal rod 
was originally installed in the NRX reactor and 
irradiated until the central portion ofthe rodhad 
received an exposure of about 3000 Mwd/tonne. 
Reference 84 describes the experimental pro- 
gram and the analytical techniques used, which 
included isotopic analysis of sections of the ir- 
radiated rod. Kushneriuk’s calculations treat the 
neutron absorption in an NRX rod, accounting for 
the effects of neutron spectrum seenat different 
cross sections of the rod andincluding the effect 
of changes in rod composition with exposure. 
The World-Average values for the nuclear pa- 
rameters—pboth the 2200 meters per second 
values and the variation with neutron energy — 
were used, with some minor revisions by 
Westcott" to include more recent valués. The 
measured initial conversion ratio (whichagrees 
with calculation) was used both in the calcula- 
tions and in the interpretation of the isotopic 
analyses. The effective cross sections most 
susceptible of accurate interpretation from the 
isotopic analyses are for the parasitic absorp- 
tions* 6, of Pu’®®, Pu’! and Pu”4!, The interpre- 
tation involves use of the fission cross sections 
6, of Pu’®® and Pu”! (and other constants); how- 
ever, for the relatively low irradiation exposures 
involved, the G, values obtainedare not sensitive 
to the 6, values (and tothe other constants) used 
within their estimated order of accuracy. For 
disks P and G of rod No. 683 (~970 and 2740 
Mwd/tonne exposure, respectively) the analyses 
give G, (239) =411+7 and 41443 barns, re- 
spectively. These compared to calculated values 
of 410 and 402 barns, respectively. Aplotis given 
of the “time-average” value of 6, (240) versus 
irradiation time since 6, (240) varies signifi- 
cantly with irradiation exposure. Kushneriuk’s 
calculated “time-average” curve agrees well 
with the data points obtained from the analyses. 
For G, (Pu’*') a value of 430 + 20 barns was ob- 
tained from disk G. No comparison with calcula- 
tions is given; however, Westcott’s tabulation” 
for the spectrum appropriate tothe NRX reactor 
(r = 0.07 where ris definedina previous section) 
also gives 430 b. Because of the over-all good 


*The symbol G designates the effective cross section 
according to the Westcott formulation (see the imme- 
. diately preceding section). 
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agreement, there is increased confidence in the 
use of the effective cross-section values for 
Pu?*?, Pu’4° and Pu’4! when proper account of 
the neutron-spectrum effects is taken. 

Another important piece of evidence relating 
to the predicted exposure values is based onthe 
results of the GLEEP pile-oscillator reactivity 
measurements for irradiated fuel samples. The 
experimental results were corrected for the 
buildup of Sm'** and for the other accumulated 
fission-product poisons. The corrected results 
then give the cross-section change, inbarns per 
initial fissile atom (bifa) due to U**> depletion 
and plutonium-isotopes production. For compar- 
ison, the reactivity variations calculated on the 
basis of effective cross-section values derived 
from Kushneriuk’s calculations for the GLEEP 
neutron spectrum are interpreted in terms of 
bifa. A large discrepancy exists between the ex- 
perimental and the calculated curves of bifa vs. 
irradiation exposure. The calculated results in- 
dicate significantly less reactivity than experi- 
ment. It is possible to force agreement by postu- 
lating “reasonable” changes in the nuclear data. 
For example, if ¢,(239)/G,(235) is increased 
by 3 per cent [with 6,(239) unchanged] and also 
v(239)/,(235) is increased by 1 per cent, the 
calculated curve comes close to the experi- 
mental curve. 

On the basis, mainly, of the foregoing evidence 
(and from other considerations discussed in the 
paper), it is concluded that the calculational 
scheme used should give a lower bound to the 
reactivity lifetime, assuming that the initial re- 
activity and resonance escape factor (and hence 
ICR) are sufficiently well established. The cal- 
culations indicate that fuel exposures greater 
than 8000 Mwd/tonne are possible in a D,O 
power reactor fueled with natural uranium oxide. 

As summarized in a previous section, the set 
of nuclear data adopted by the British*!.® gave 
good agreement”? of calculated and measured re- 
activity changes with irradiation. The use ofthe 
World-Consistent set of nuclear data (BNL-325 
and its Supplement 1 plus Westcott’s revisions) "4 
appears to underestimate the contributionto re- 
activity due to the buildup of the plutonium iso- 
topes with irradiation exposure; this contribution 
to reactivity for a given initial conversion ratio 
depends on the relative magnitudes of the nuclear 
parameters for the plutonium isotopes and aaa 
However, rather small changes in the nuclear 
data are sufficient to give good agreement. This 
behavior emphasizes the delicacy of the balance 





between the neutron credits and debits which 
determine the long-term reactivity change, and 
the importance of accuracy in the nuclear con- 
stants if those changes are to be predicted by 
calculation. 
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THERMAL CONSIDERATIONS 





Heat Transfer and Thermal 
Cycles for Gas-cooled Systems 


The gas-cooled power reactor, which fora num- 
ber of years saw little development inthe United 
States as a central-station type, is now becoming 
an important part of the U.S. program. The heat- 
transfer characteristics of gaseous coolants are 
quite different from those of liquids, and the 
interrelationships among the thermal variables 
to determine an optimum plant designare rather 
different from those which apply for liquid- 
cooled reactors. Consequently, studies which 
treat the relationships of the thermal variables 
to each other and to other aspects of the plant 
design are particularly pertinent at this stage of 
the gas reactor program inthe United States. The 
papers reviewed below reflect mainly the results 
of the British investigations in these areas. 


Probably the outstanding features of gas- 
cooled systems as distinguished from those 
cooled by liquids are the lower jheat-transfer 
coefficients, the large coolant volumes required, 
and the rather large amounts of power required 
to circulate the gas coolant. The power consumed 
in circulating the coolant ina gas-cooled reactor 
plant is typically of the order of 8 per cent of 
the plant’s gross electric output. It is therefore 
very important to minimize the pressure drop in 
the gas-coolant circuit. 


A large part of the power consumption con- 
nected with coolant circulation occurs in the 
steam generators of a gas-cooled power reactor 
system. Reference 1 points out that actually two 
major consumption items occur in the steam 
generator. The first is the gas side pumping 
power associated with the gas pressure drop 
through the steam generator heat-transfer sur- 
face. The second is the power loss implied by 
the existence of finite temperature differences 
between the gas and water sides. This makes the 
transfer of heat irreversible and therefore 
causes a loss of theoretically available power. 
The two power losses can be summed to obtain 


the total power loss incurred by a steam gener- 
ator. 

The power loss due to the temperature differ- 
ence can be increased or decreased by changing 
the amount of heat-transfer surface. The change 
in heat-transfer surface can be obtained by 
changing either the heat-exchanger face areaor 
its length. In either case. the gas pressure drop, 
and therefore the power loss due to pumping 
power will be affected. Reference 1 points out 
that there is an optimum balance between the 
heat-transfer performance, as reflected in the 
gas-steam temperature difference, and the pres- 
sure drop for the steam generators of a gas- 
cooled reactor. 

The reference also presents a method of evalu- 
ating and correlating these interrelated parame- 
ters for any given type of heat-transfer surface 
configuration. The graphically illustrated re- 
sults presented are based on the circumferen- 
tially finned steel Integron tube used in the 
British Bradwell power station. They are not in- 
tended to cover heat-transfer tube dimensions 
or arrangements differing greatly from that 
particular type. Qualitatively similar results 
could, however, be obtained for other configura- 
tions. 

The graphical results presented in the refer- 
ence include only the thermodynamic considera- 
tions. No consideration is given to the capital 
investment costs involved, but importance of this 
factor is acknowledged. The effect of taking the 
cost of heat-transfer surface into account would 
be to shift the optimum parameters toward less 
heat-transfer surface. 

A more complex problem in gas-cooled reac- 
tor design is the optimization of the reactor 
coolant gas conditions and the steam cycle con- 
ditions for lowest power cost. 

Reference 2 presents a study that evaluates 
the thermodynamic optimum steam cycles and 
their efficiencies for gas-cooled reactors as the 
gas outlet temperature is increased. The range 
of gas outlet temperatures considered is from 
the 662°F of Calder Hall to 1110°F. The types of 
steam cycles considered are single pressure, 
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dual pressure, and super critical, each withand 
without reheat. The analysis is based on reac- 
tors of the Calder Hall type. 

Table III-1 shows steam cycles and operating 
conditions that were found to be optimum as the 
attainable outlet coolant temperature was as- 
sumed to increase. It is to be notedthat the sin- 
gle pressure cycle is not indicated to be optimum 
at any point over the temperature range consid- 
ered. The dual-pressure cycle showed the 
highest gains over the single-pressure cycle 
when the coolant temperature leaving the steam 
raising plant was low. The relative positions of 
the two cycles would no doubt be changed if the 
increased cost of equipment for the dual- 
pressure cycle, especially that of the heat ex- 
changers. were taken into account. 


inlet temperature can be increased to the point 
that feed-heating can be used to great advantage. 
Much of the potential improvement in cycle effi- 
ciency with increased coolant outlet temperature 
stems from the fact that the coolant inlet tem- 
perature can be raised as well. The restriction 
on inlet temperature for a given outlet tempera- 
ture arises from pumping-power considerations 
and from the circumstance that some specified 
amount of power must be removed from the re- 
actor. In the analysis of reference 2 the pumping 
power is held constant at 10 per cent of the net 
power output of the plant, and the reactor power 
density is evidently held constant, at an unspeci- 
fied value. 

The report apparently evaluates cycles en- 
tirely on a thermodynamic basis. Since no men- 


Table II-1 ADVANCED REACTOR STEAM CYCLE DATA 





Coolant outlet temperature, °C 350 
Coolant outlet temperature, (°F) (660) 
Coolant inlet temperature,* °C 177 
Coolant inlet temperature, (°F) (350) 
Type of steam cycle Dual pressure 
H. P, steam pressure, psi 600 
H. P. steam temperature, °C 640 
H. P, steam temperature, (°F) (337) 
L. P. steam pressure, psi 150 
L. P. steam temperature, °F 620 
L. P. steam temperature, [°C) (327) 


Reheat pressure, psi 
Reheat temperature, °F 
Reheat temperature, (°C) 


Proportion of h.p. steam, % 60 
Feed water temperature, °F 200 
Feed water temperature, (°C) (93) 
Turbine cycle efficiency, % 28.2 
Turbine exhaust moisture, % 13.5 
Station over-all efficiency, % 25 
Efficiency with single- 

pressure steam cycle, % 23.0 
Corresponding turbine 

exhaust moisture, % 11.0 








450 500 550 600 
(840) (930) (1,020) (1,110) 
260 294 320 343 
(500) (560) (610) (650) 
—_—_—_—_— —. v 4 
Dual pressure Super -critical 
with reheat with reheat 
2,000 2,800 3,500 5,500 
820 910 1,000 1,090 
(438) (487) (537) (590) 
500 860 
820 910 
(438) (487) 
1,200 1,800 
1,000 1,090 
(537) (590) 
67 79 100 100 
455 515 520 580 
(235) (268) (272) (305) 
38.0 41.0 43.5 45.0 
10.5 11.5 11.5 12.0 
34 37 39 40.5 
31.5 34.5 36.5 37.5 
14.5 t t Tt 





*Coolant temperature at inlet to circulators. 
TExcessive moisture content. 


A large increase in the attainable steam pres- 
sure and a correspondingly attractive improve- 
ment in the over-all thermal efficiency are seen 
to result from relatively modest advances in re- 
actor coolant outlet temperature. When coolant 
outlet temperatures are limited to about 750°F, 
the effective use of feed-heating is not possible 
since it would be done at the expense of steam 
pressure. As the coolant outlet temperature is 
increased above this value, however, the coolant 


tion is made of capital investment costs, the 
method of selecting the optimum cycle is not 
clear. It appears that the optimum cycles se- 
lected are very near those which give maximum 
over-all plant efficiency with an acceptable tur- 
bine exhaust moisture. The report does point out 
that in the final selection of optimum cycle con- 
ditions, practical and economic considerations 
must supervene and that economic considera- 
tions lead to selection of coolant return tem- 
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peratures (and hence efficiencies) lower than the 
indicated thermodynamic optimum. As an ex- 
ample of this, the thermodynamic optimum indi- 
cated by the results of the report would .corre- 
spond to a coolant temperature at the exitof the 
steam raising plant of 370to 400°F for the 662°F 
reactor exit temperature of Calder Hall. In the 
feasibility study for Calder Hall, however, the 
coolant temperature leaving the steam raising 
plant was selected at 350°F. 


In the foregoing discussion no details of the 
heat transfer within the reactor were considered. 
The internal heat transfer is of great importance, 
for it affects the total power attainable from a 
given reactor, the pumping power, and the attain- 
able outlet temperature. Extended heat-transfer 
surfaces are oftenused in gas-cooled reactors to 
obtain higher heat fluxes from the fuelfora given 
temperature difference between fuel and coolant. 
A discussion of earlier data on the Calder Hall 
transverse fins was given in a previous issue of 
this review.’ The British have also performed 
experimental work onother configurations of ex- 
tended surfaces. Reference 4 presents test re- 
sults on a configuration consisting of multistart 
helical fins intersected by four longitudinal 
splitters. Correlations of heat transfer and fric- 
tion factor were made relative to the number of 
starts and the lead of the helix over a range of 
Reynolds number. 

The dimensions of the finned tube configuration 
treated in reference 4, which were constant for 
all cases considered, were 


Outer diameter of fins, in. 2.0 
Diameter of fin root, in. 1.25 
Number of splitters 4 
Thickness of splitters, in. 0.064 
Radial clearance between outside 

of splitters and channel wall, in. 0.050 


Finning material Magnox A-12 
The range of variables considered was 


Inside diameter of flow 


channel, in. 3.5 to 4.125 
Number of helix starts 30 to 48 
Helix leads, in, 12 to 40 
Metal to gas temperature 

difference, °C 50 to 90 


Reynolds number 3 x 104 to 7 x 104 


Atmospheric air was used as the coolant gas for 
all tests. Temperature dependent variables were 





determined by using average thermocouple read- 
ings. 

The test results obtained show that the follow- 
ing equation predicts the Stanton number and 
therefore the heat-transfer properties to +12 
per cent within the 95 per cent confidence limits. 


8800 2500 
=] 6. Pee See 
St = 10 ( N + L 44) 


"= Se (L-22) 
16 


where N= number of fin starts 
L = lead of helix, in. 
kg and k m = thermal conductivity of gas and 
metal, respectively 


The Stanton number was found to vary about 
+10 per cent ina definite sinusoidal pattern along 
the length of the test section. The friction factor, 
as calculated from pressure-drop measure- 
ments, includes the effect of mixing of the gas by 
the fins. The following equation was found to 
correlate the data within +20 per cent. 


f = Re~?- 151 ne + 0.0126) 


Gas velocity profiles were taken over the flow 
cross sections. These showed local velocities of 
about 80 per cent of the mean along the tube and 
fin surfaces and as high as 155 per cent of the 
mean in the approximate centers of the flow 
sections. With anincrease in flow channel diame- 
ter, and constant tube and fin dimensions, it 
might be expected that the section bounded by the 
splitters, the fin tips, and the channel wall would 
rob the gas flow between the fins and thereby 
significantly reduce the effective heat-transfer 
coefficient. Comparison of results for different 
diameters showed, however, that such an effect 
did not occur over the interesting range of di- 
ameters. 

The data reported in reference 4 for two dif- 
ferent flow channel dimensions are as follows: 


€ eS 2S ae aes 
3.5 0.695 30 18 515 1.0 1.47 
4.125 0.932 30 18 509 1.47 1.0 


where C = flow channel diameter, in. 
D, = equivalent hydraulic diameter, in. 








26 GENERAL RESEARCH AND DEVELOPMENT 


N = number of fin starts 
L = length of helix lead of fins, in. 

St, = Stanton number for: (Reynolds num- 
ber)(ratio of gas conductivity/fin 
conductivity) = 8 

Af¢= relative flow area 

V = relative velocity 


The theoretical effect of diameter onthe heat- 
transfer coefficient, h, assuming no by-passing 
of the flow, should be given by the equation: 


Nusselt number = 
(constant) (Reynolds)®:® (Prandt1)°-4 


hy. 125 = (Ys.ne)"" x ( Ds 5 - = 0.69 
hs.5 V3.5 Dg 125 





{f the relative values of h based upon the meas- 
ured values of Stanton number (St, above) are 
compared: 





Stanton number =h te ") 


here C, the specific heat, and p, the density of 
the gas, remain constant for the two cases. Then 


Thus, for this geometry and flow condition, the 
calculated ratio of heat-transfer coefficients, 
using average conditions and the conventional 
heat-transfer equation, is very close to the 
measured ratio. It therefore may be said that no 
robbing of flow effect was observed. The obvious 
conclusion is that good mixing is provided by the 
splitters and the helix. 

In the previous discussion of the Calder Hall 
fins in this Review*® plots were made of Nusselt 
number vs. Reynolds number and of (Stanton 
number) (Prandtl number)°-°/('/, friction factor) 
vs. Reynolds number; the latter curve gives the 
heat-transfer/friction ratio, and is anindication 
of the price paid in friction pressure drop (and 
therefore pumping power) for the improved heat 
transfer provided by the fins. The results of 
reference 4 can be compared inthis respect with 
the Calder Hall results. If the value of the ratio 
of k (gas)/k (magnox fins) is takenas about 2.6 x 
10~*, then the St, measured for the extended 
surface configuration tested in reference 4cor- 
responds to a Reynolds number of 3.09 = 10‘, and 


St;, corresponds to Re = 6.18 x 10‘. The corre- 
sponding Nusselt numbers are 110 and 170, re- 
spect:vely. When these values are plotted on 
the previously mentioned graphs* of reference 3, 
it may be seen that the extended surface configu- 
ration of reference 4 gives a performance about 
equal to that of the “roughened tube.”’ 

When heat-transfer tests are made on systems 
ccntaining extended surface, the direct resultis 
an “effective’’ heat-transfer surface coeffi- 
cient. Often the desired result is a “local’’ or 
true surface coefficient. Reference 5 presentsa 
method of directly evaluating the “local’’ coeffi- 
cient without a trial-and-error solution that is 
required by the more straightforward procedure. 
The method presented is an approximate solution 
that may introduce an error of about 2 per cent. 
In what is probably the more typical situation of 
performing design calculations, rather than ana- 
lyzing test data, the problem discussed by this 
paper does not occur. 

Reference 6 derives the appropriate mathe- 
matical equations for use in analogue computer 
studies of the transient characteristics of a gas- 
cooled graphite-moderated reactor. 


Properties of Gases 


Problems involving the properties of gases 
occur in reactor technology for widely different 
applications such as gas cooling, energy absorp- 
tion in containment gases, gas-bonding of fuel and 
jackets, and hypothetical reactor explosions. 
Fairly accurate (<5 per cent error average) 
calculations for thermal conductivities of gas 
mixtures can be made, with quite limited pure- 
gas input data, by using formulas derived in ref- 
erences 7 and 8. These formulas are correlated 
with experimental data. 

Reference 7 is a detailed thesis with new, 
careful experiments to determine thermal con- 
ductivity and viscosity for pure and mixed gases 
at temperatures up to 540°C. The empirical for- 
mula offered by reference 7 fits experimental 
data a bit more accurately than the quite similar 
expression of reference 8 which has the advan- 
tage of having one less input parameter (that is 
difficult to determine) and further, is derived 
from rigorous theory with certain approxima- 
tions to reduce the complexity of calculations. 
Obvious advantages occur for use of this formu- 


*Figs. 4 and 5, pp. 9 and 10, of reference 3. 
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lation when high-temperature problems are con- 
sidered for gas mixtures in connection with deto- 
nations and shock waves. 

When mixtures of fission gases (xenon and 
krypton) with helium are considered, the large 
differences in atomic weight cause marked devi- 
ations from the simple mixing law. A helium 
“bond gas,’’ used to improve heat transfer from 
a fuel element like UO, to its metallic jacket, 
will have a slow dilution from krypton and xenon 
with exposure since some fission gases escape 
the solid fuel and mix with the helium. The 
relatively low thermal conductivities of xenon 
and krypton seriously impair the conductivity of 
the gas bond, and itis important to allow for this 
effect in the fuel element design. The formulas of 
references 7 and 8 are applicable and may be 
found useful for this case. 

Reference 9 concerns essentially the same 
problems faced by references 7 and 8 but for 
considerably higher temperatures, 1000°K to 
15,000°K, which are beyond the present realm of 
laboratory tests. The basic theory of molecular 
interactions is defined using experimental re- 
sults from molecular-beam-scattering experi- 
ments. The results offer anapproach only, since 
effects of excitation, dissociation, and ionization 
are not considered. 
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REACTOR CONTROL: BURNABLE POISONS 





The inclusion of a burnable poison in reactors 
has been considered for some time as a means 
of extending the reactivity lifetime and de- 
creasing the number of conventional control 
rods required. Burnable poison was first used 
in the BORAX-III reactor,’ in 1955. In reactors 
having conversion ratios less than unity, the 
burn-up of fissionable material gradually de- 
creases the available reactivity in the reactor 
during the life of the fuel loading, until even- 
tually the available reactivity is not sufficient 
for operation of the reactor. The length of time 
required to reach this point may be adjusted by 
the amount of initial excess reactivity incor- 
porated into the fuel loading. However, in the 
early stages of the reactivity lifetime, this ex- 
cess reactivity must be controlled by the pres- 
ence of control rods or other control means. 
A burnable poison may be used as one of these 
means of control. It consists of some neutron- 
absorbing isotope which, upon the absorption of 
neutrons, is converted to an isotope of low ab- 
sorption. For this purpose materials are chosen 
of such cross sections that the rate of burn-up 
of the absorbing isotope approximately compen- 
sates for the rate of reactivity loss by fuel 
burn-up. Two advantages which result from the 
use of a burnable poison for the added control 
are (1) a smaller flux distortion than in the 
case of control rods and (2) a saving in the 
number of control rods required and a saving 
in space both in the control rod drive area and 
in the reactor core. 

The usual application of burnable poisons is 
in highly enriched reactors containing negli- 
gible quantities of fertile materials. For ther- 
mal reactors of this type, if the fission-product 
poisoning and the change in neutron leakage with 
diffusion length are neglected, the relation be- 
tween reactivity change and fuel burn-up is 
approximately 


en. (1 - 5) 


k est M q 


where dM/M is the fractional change in the 
mass of fissionable isotope in the reactor, k is 
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the infinite multiplication constant, and 7 is the 
number of neutrons produced per neutron ab- 
sorbed in the fissionable isotope. The range of 
(1 — k/n) for practical reactors is perhaps from 
about u.1 (k = 1.86) to about 0.5 (k = 1.04), and 
in many fully enriched reactors it is between 
0.3 and 0.4. It is clear that for most such re- 
actors the change in reactivity withtuel burn-up 
is rapid, and if, for example, it is desired to 
burn 50 per cent of the fuel in the core, the 
corresponding reactivity change will be very 
large unless some scheme such as burnable 
poisons is used. 

Ideally there is a certain effective cross 
section for the burnable poison which will cause 
it to burn out at a rate which very nearly com- 
pensates for the reactivity loss due to fuel 
burn-up and fission-product poisoning. Further- 
more, it is desirable that the burnable poison 
be completely consumed so as to minimize the 
total fuel inventory required. Thai this ideal 
behavior cannot be expected from a single uni- 
formly distributed burnable poison is apparent 
from rather simple considerations. The equation 
for Kerf of a uniform, thermal, nonconverting 
reactor can be written as 


k eff = - : 
eff 1 + M’B? 1 + (5,/2:5) 





where 2; is the macroscopic absorption cross 
section of the fissionable isotope and £, is the 
macroscopic absorption cross section of all 
other materials in the reactor (including the 
burnable poison). If again fission-product poi- 
soning and the change in neutron leakage (i.e., 
the change in M’B’) are neglected, then keg 
would remain constant with fuel burn-up pro- 
vided the ratio Z,/2.; remained constant. If the 
burnable poison were the only absorber present 
other than the fissionable isotope, the ratio 
would remain constant if (and only if) the 
microscopic absorption cross section of the 
burnable poison were equal to that of the fis- 
sionable isotope. However, the poison would not 
4e completely burned out until the fissionable 
isotope was also completely burned out. In 
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practice, Zr is made up of other absorptions 
besides that of the burnable poison. For this 
case, if one can choose any desired micro- 
scopic cross section for the burnable poison, 
one can make the ratio £,/Z,, equal to its 
initial value at one chosen degree of fuel 
burn-up, but at no more than one. To do this, 
the fractional rate of poison burn-up must be 
chosen to be higher than that of the fuel by an 
appropriate amount (that is to say, the micro- 
scopic absorption cross section of the poison 
must be larger than that of the fuel). Obviously 
the initial rate of burnout of the poison will then 
be too high. Thus the characteristic behavicr of 
reactivity when a single, uniformly distributed 
burnable poison is used is an initial increase 
with fuel burn-up, followed by a decrease which 
eventually returns the reactivity to its initial 
value, followed by a rather rapid reactivity loss 
with further burn-up. 

Numerous and more complicated schemes 
have been devised to circumvent the difficulties 
indicated above. The analysis of these schemes 
is complex, and the flux perturbations caused by 
control rods and by the burnable poison itself, 
as well as the changing neutron spectrum, 
complicate further the calculations involved in 
determining the optimum poison concentration 
and configuration. Some of the burnable poison 
schemes and the methods of Calculation are 
given in references 2 and 3. 

The most commonly considered material for 
a burnable poison is boron, although dysprosium 
has also been studied. In epithermal reactors, 
cadmium shows some promise, although the re- 
sults are probably quite sensitive to the reactor 
spectrum. Both boron and dysprosium have 
thermal absorption cross sections considerably 
greater than that of U*®*. The effective cross 
section may be adjusted, however, by taking 
advantage of the phenomenon of self-shielding in 
lumped poison rods. With heavy self-shielding, 
the effectiveness of the poison does not change 
with burn-up during the early period of ir- 
radiation, and during this period the reactivity 
of the reactor decreases as the U**® is burned. 
Later, after enough of the poison is consumed 
to reduce the self-shielding, the burnable poison 
is consumed faster than the fuel, and the re- 
activity rises to a maximum before again 


dropping off with further irradiation. Radkowsky’ 
has proposed that the minimum mis-match be- 
tween poison burn-up and fuel burn-up may be 
obtained by using burnable poison elements of 
two different self-shielding factors, adjusted 
so that the reactivity peak of the lightly shielded 
poison element coincides with the reactivity dip 
of the heavily shielded poison element. 

Little has been written on the applicability of 
burnable poisons to converter reactors. It would 
be expected that in most cases they would not 
be attractive, since their use implies the con- 
sumption of neutrons which might better be 
used for conversion. This reduction of neutron 
economy is no greater however than would re- 
sult from the use of control rods for the same 
purpose, provided the burnable poison is vir- 
tually all consumed during the life of the fuel 
charge. In the case of the U?** converter the 
behavior of the single, uniformly distributed 
burnable poison appears unattractive because 
the reactivity peak provided by the burnable 
poison would seem to coincide with the plutonium 
reactivity peak. In thorium converters this con- 
sideration does not hold, and indeed, it is pos- 
sible that a high cross-section burnable poison 
might serve to level out the initial reactivity 
dip due to the hold-up of potential U’*™ as 
protactinium. It has been stated that burnable 
poison is used in the Consolidated Edison tho- 
rium converter reactor.‘ 
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SHIELDING 





Neutron-flux Distribution 
and Attenuation 


In a paper presented at the second Geneva con- 
ference,' Blizard outlines reactor shielding and 
discusses basic data requirements, facilities for 
experimental research, attenuation theory, and 
recent developments in shielding materials. 
There is a growing need for better neutron 
cross-section data and for improved methods to 
determine neutron flux distributions in shields. 

Most of the experiments to determine the 
widely used fast neutron removal cross sections 
have been performed at the ORNL Lid Tank 
Shielding Facility (LTSF). The LTSF measure- 
ments require large (142 by 167 cm) slabs of 
the material under test; for some materials 
these quantities may be expensive. Reference 2 
describes an experiment utilizing relatively 
small samples of material. The experiment was 
performed at the Brookhaven Shielding Facility 
with an air-filled aluminum collimator placed 
above the converter plate and surrounded by 
water to provide a fairly small diameter source 
of essentially monodirectional neutrons. The 
results of this experiment yield “collimated” 
removal cross sections for hydrogen, lithium, 
beryllium, boron, carbon, oxygen, fluorine, 
sodium, aluminum, sulfur, chlorine, titanium, 
iron, nickel, copper, cadmium, tungsten, lead, 
and bismuth. These cross sections are sig- 
nificantly larger than the conventional removal 
cross sections obtained with a sample close to 
a distributed isotropic source but can be used 
for an indication of the relative effectiveness of 
different shield materials. The monodirectional 
nature of the source used in the experiment 
prevents use of the results in the usual way in 
shield design calculations. 

Trubey and Chapman’ have used the ORNL 
LTSF to obtain removal cross sections for car- 
bon and oxygen when the material is distributed 
rather than concentrated adjacent to the source 
plate. The carbon cross section was determined 
from a comparison of thermal neutron flux meas- 
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urements in a sugar solution with measurements 
in plain water. Since the sugar (C,.H,.0;;) 
solution contained almost as much hydrogen 
and oxygen as plain water and contained them 
in the same ratio, the effect of the water 
could be separated and a distributed cross 
section for carbon obtained. Then by com- 
paring the measurements in oil and water 
and using the carbon cross section found as 
described above, the oxygen cross section was 
determined. Though differences appear between 
these “distributed” cross sections andthe cross 
sections determined by the usual method, the 
indicated differences are well within the ex- 
perimental uncertainties involved and lead to 
the conclusion that any actual differences are 
probably small. 

Bourgeois et al.‘ present removal cross sec- 
tions obtained experimentally and theoretically. 
The majority of the values given are the results 
of experiments performedat the French NAIADE 
facility and agree fairly well with previously 
reported cross sections. The theoretical cal- 
culations treated energy-dependent fast neutron 
removal cross sections for iron, lead, and 
water. 

A British paper’ presents a method for pre- 
dicting the spatial distribution of thermal and 
epithermal neutron fluxes in hydrogen-metal 
reactor shields. Such calculations are particu- 
larly applicable to the prediction of heat pro- 
duction in thermal shields, pressure vessels, 
and other materials relatively near the reactor 
core. The three-group model described provides 
a value of the thermal flux accurate to better 
than 20 per cent. Empirical parameters are 
required for application of the method to com- 
plex multilayered systems. The experiments 
were performed at the Harwell shielding re- 
actor, LIDO, which consists of a core of en- 
riched U-Al alloy suspended in a water-filled 
concrete tank with aluminum-faced cavities at 
three positions around the tank walls. Large- 
scale shielding experiments may be constructed 
outside a cavity and the core may thenbe moved 
to within a few centimeters of the inner face of 
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the cavity. Iron was chosen as the metal in the 
hydrogen-metal shield, and after empirically 
deriving the slowing-down lengths in the shield 
materials, applying these lengths to the three- 
group model, and comparing the semianalytic 
results with further experimental measure- 
ments, 20 per cent or better accuracy for the 
thermal flux distribution was found to prevail 
in shield systems containing from 2.4 to 37.4 
cm of iron. 

Duncan and Whittum® have developed a rela- 
tively simple method of determining fast and 
thermal neutron flux distributions at deep neu- 
tron penetrations in hydrogenous shields. The 
method uses fast neutron removal theory in 
conjunction with a two-group diffusion model 
for calculating the thermal neutron flux distri- 
bution. It is accurate only in hydrogenous shield 
regions several diffusion lengths from non- 
hydrogenous boundaries. The authors note that 
since the method overestimates the thermal flux 
in the nonhydrogenous portions of a shield and 
is simple in comparison with other methods, it 
may find use in preliminary design studies. 


Concretes for Reactor Shielding 


The low unit cost of concrete has resulted in 
its use as shielding for nearly all power re- 
actors where shield thickness and weight are 
not critical. The “standard” methods have gen- 
erally been employed to calculate required 
thicknesses of concrete; i.e., the moments 
method for gamma-ray attenuation and the re- 
moval cross section concept for fast-neutron 
attenuation. The removal cross-section concept 
was originally developed for hydrogen-rich 
shields. It assumes a scatter collision to be 
equivalent to an absorption. It would be ex- 
pected to apply only to shields which have a 
relatively short slowing-down length for first- 
collided neutrons. Since the validity of the con- 
cept thus depends primarily on the slowing-down 
capability of the shield-contained moderator, 
its applicability to concrete shields without added 
hydrogen has been questioned. However, con- 
crete contains a considerable amount of oxygen 
which by elastic scattering augment the hydrogen 
moderation. 


To determine the actual hydrogen require- 
ments for concrete shielding and to compare 
measured gamma-ray attenuation to moments- 
method calculations, Blizard and Miller’ per- 


formed a series of tests on ordinary concrete 
at the ORNL LTSF. The rather thorough tests 
show that with no added hydrogen, ordinary 
concrete with a 7 wt.% water content attenuates 
the intermediate-energy neutrons considerably 
more strongly than it does the very fast neu- 
trons, thereby indicating adequate moderation. 
The conclusion is that in concretes containing 
at least 7 wt.% water, the shield attenuation rate 
is characteristic of the very fast neutrons and 
is adequately described by the removal cross- 
section method. The measured gamma-ray at- 
tenuation was in good agreement with calculated 
attenuation. 

Although the test procedure is not described, 
reference 4 contains experimental values for 
removal cross sections of concrete obtained 
in the French NAIADE device, a facility very 
similar to the ORNL LTSF. Ordinary and heavy 
concretes were tested, and the cross sections 
obtained Zp, = 0.085 cm-! to 0.125 cm~! are 
generally in good agreement with those reported 
in references 7 and 8. The aggregates utilized 
in the French concretes result in mixes chemi- 
cally similar to those used in the United 
States." 8.9 

The measured attenuation inthe ORNL Graph- 
ite Reactor shield is described in detail by 
Blosser et al.,® and cross plots of the data 
should furnish interesting information concern- 
ing radiation streaming through a cylindrica 
void in a concrete shield. 

The articles reprinted in reference 9 deal 
principally with selection and proportioning of 
the constituents of concretes, but also yield 
the following information: (1) barytes concrete 
mixed and placed by conventional methods has 
about twice the compressive strength of pre- 
packed barytes concrete; (2) conventional ba- 
rytes concrete compares favorably with ordinary 
concrete when subjected to severe thermal cycle, 
but prepacked barytes concrete is inferior; and 
(3) conventional barytes concrete and prepacked 
barytes concrete exhibit almost twice the ther- 
mal coefficient of expansion of ordinary con- 
crete and have relatively very low values of 
specific heat, thermal conductivity, and dif- 
fusivity. Reference 8 indicates that ordinary 
and barytes-haydite concretes suffer about the 
same percentage loss in compressive strength 
when subjected to heat. 

Clark’® has published a compilation of data 
obtained by a literature survey on the subject 
of radiation damage to concrete. The only ap- 
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parent deleterious effect resulting from irra- 
diation is the production of gas by boron- 
containing cements. For the majority of 
reactors, it is unlikely that there would be any 
appreciable change in the properties of the 
shield due to radiation damage during the 
operational lifetime. Afar more serious problem 
is overstressing as a result of nuclear and 
external heating. 
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FUEL ELEMENTS 





Electron-beam 
Vacuum -welding Process 


Because of the potential economic and technical 
advantages of electron-beam welding, this weld- 
ing technique is being investigated for making 
fuel-element closures and welding reactor struc- 
tural materials. At present, the most commonly 
used welding method for these applications is 
the inert-gas shielded tungsten-arc technique. A 
disadvantage of this method is the possible weld 
contamination from either the electrode mate- 
rials or impurities in the inert-gas shielding 
atmos: ‘ere. Difficulty is also encountered inthe 
joining of components of dissimilar thickness. 

The vacuum-welding technique employs a 
highly concentrated beam of rapidly moving 
electrons which locally fuses the parts to be 
joined. The electrons are emitted by a heated 
tungsten filament, accelerated through a high 
electrostatic potential difference, and focu3ed 
electrostatically to form a “spot’”’ from 1; to % 
in. in diameter at the junction to be welded. The 
size and penetration of the molten zone are func- 
tions of the electric power fed into the electron 
beam and of the degree of focusing. Welding is 
done at a pressure of 107° mm Hg or less. At 
higher pressures, ionization canoccur, resulting 
in arcing and consequent loss of focusing ability. 

The electron-beam welding process is saidto 
possess a number of significant advantages over 
the customary shielded arc technique. First, the 
fusion welding of widely differing thicknesses of 
materials can be done easily. This is accom- 
plished by adjusting the focus of the beam to ob- 
tain the desired “spot’’ diameter and the accel- 
erating potential to give the required penetration. 

Hanford’ has been developing the electron- 
beam vacuum-welding process to overcome 
some of the limitations of the inert-gas shielded 
process. The method was applied to the welding 
of end caps on Zircaloy-2-, aluminum-, and 


stainless-steel-clad fuel elements and to other 
component materials such as copper, gold, 
nickel, silver, platinum, and uranium. 

Hanford! has reported that a power input of 
20 watts is sufficient to weld foil, whereas up to 
1000 watts is used to fuse satisfactorily '4-in.- 
thick materials. It was also claimed that very 
low contamination of the material in the weld 
zone was achieved. Since no electrode is used, 
contamination of the weld by electrode materials 
is impossible. Similarly, the high vacuum re- 
quired by the process reduces impurities in 
the welding atmosphere to a negligible quantity. 
The result is a product having more ductility 
and greater corrosion resistance than can be 
obtained with the more conventional welding 
processes. Corrosion tests of electron-beam 
vacuum-welded Zircaloy-2 fuel cans in 686°F 
deionized water indicated that corrosion resist- 
ance of the weld and heat-affected zone was as 
good as or better than the parent metal. 

Another advantage mentioned in the reference 
report is the improved mixing of the molten ma- 
terial, to yield a more uniform metallurgical 
structure in the fused zone. The electron-beam 
welded specimens are characterized by a homo- 
geneous acicular grain structure with a com- 
plete absence of gas pockets or “pipes’’ of 
foreign materials. Tests are reported! to indi- 
cate that the properties of the fused and heat- 
affected zones are superior to those of the parent 
material. 

The electron-beam welding process has been 
used successfully by the French.? Aluminum- 
clad fuel elements, for the Saclay EL-3 reactor, 
had 2-mm-thick end caps welded to 1-mm-thick 
cans by this method. Stainless steel, zirconium, 
beryllium, aluminum, molybdenum, uranium, 
and magnesium were also welded using this 
technique. Some difficulty was encountered by 
the French with the welding of magnesium be- 
cause the magnesium vapor formed during the 
fusing led to ionization and arcing inthe vacuum 
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chamber. This difficulty was overcome by using 
a pulsating source of electrons and proper con- 
trol of the motion of the piece being welded. A 
large size electron-beam vacuum-welding ma- 
chine has been constructed in France to weld 
large zirconium assemblies. This machine was 
used to weld the tank of “Proserpine,’’ a pluto- 
nium salt-fueled homogeneous reactor.’ 

Bettis Plant is also constructing an electron- 
beam welding apparatus to be used to weld 
Zircaloy-2 components for the second core of 
the Shippingport Pressurized Water Reactor.’ 

Operational costs (not including labor) are 
said to be substantially lower for the electron- 
beam process than for the shielded-arc process. 
Because the electron-beam process must be 
virtually completely automatic, the initial cost 
of equipment may be high compared to that for 
the tungsten-arc method. However, for the mass 
production of fuel elements a high degree of 
mechanization is desirable, in any case, andthe 
potentially good controllability of the electron- 
beam process along with its other advantages 
may make it the economically preferable one. 


Electrocladding of Fuel Elements 


Aluminum-clad uranium-metal fuel elements 
have been used extensively in reactors. Most of 
these elements have been made up of uranium- 
metal fuel, AlSi bonded to aluminum cans‘ or 
have been U-Al alloys or dispersions roll clad 
with aluminum. Atomics International? has re- 
cently reported on their development of an elec- 
trocladding technique of canning fuel elements. 
The method involves the use of a continuously 
regenerative, nonaqueous, aluminum electro- 
deposition process. Anintermediate nickel layer 
on the uranium metal is used to promote bonding 
without formation of a brittle intermetallic bond 
zone. It was reported by Atomics International 
that dense, adherent claddings of aluminum up 
to 25-mil thickness were produced on nickel- 
plated cores by electrodeposition; metallurgical 
bonding of the electrocladdings was obtained by 
a thermopneumatic process. 

The production of an integral clad by this 
method appears to be a distinct advantage over 
methods requiring end-cap welding, particularly 
where complicated fuel-element shapes are re- 
quired such as tubular or twisted ribbon ele- 
ments. 


Because the scope of work on this project was 
limited, no economic comparison of the fabrica- 
tion methods can be made. 


Lead Bonding of Fuel Elements 


In most proposed power reactors employing 
oxide fuel, the oxide is bonded to the fuel canning 
material by filling the annular gap with helium. 
Close tolerances or swaging techniques are 
therefore required to obtaina very small gap and 
thereby maintain a small temperature drop 
across the helium gap. In addition to this, re- 
leased fission product gases mix with the helium 
and cause a drop in the thermal conductivity of 
the mixture, thereby resulting in a higher tem- 
perature drop andan increase inthe central tem- 
perature of the fuel material for the same oper- 
ating heat flux. 

A possible substitute for the helium bond is a 
lead (Pb) bond which has a relatively high 
thermal conductivity. This has been done in the 
aluminum canned ThO,-UO,, BORAX-IV ele- 
ments. However, most of the power reactors 
being designed use Zircaloy-2 or stainless steel 
for the oxide fuel can andthe compatibility of Pb 
with these materials is questionable. 

Argonne National Laboratory® has recently 
reported preliminary data on the reactions of Pb 
with ThO,-UO,, stainless steel, and Zircaloy-2 at 
temperatures of 600 and 800°C. Metallographic 
examination of the reaction zones of the speci- 
mens held at these temperatures for 60 days 
showed: (1) no reaction between Pb and ThO,- 
UO, or between Pb and stainless steel at tem- 
peratures of 600 and 800°C; and (2) no reaction 
between Pb and Zircaloy-2 at 600°C and some 
reaction between Pb and Zircaloy-2 at 800°C. 

It would appear that, in water type reactors 
where the bonding layer would probably be less 
than 350°C, a Pb or possibly a Pballoy bond for 
oxide type fuels may be feasible. 


References 


1. W. L. Wyman and W. I. Steinkamp, Electron Beam 
Vacuum Welding— Development of. Process, HW- 
55667, Apr. 8, 1958. (Unclassified AEC report.) 

2. J. Briola, Automatic Welding of Fuel Elements, 
Second United Nations International Conference on 
the Peaceful Uses of Atomic Energy, Geneva, 





FUEL ELEMENTS 35 


September 1958, A/CONF.15/P/1161. (Available 
from the Office of Technical Services, Dept. of 
Commerce, Washington 25, D. C.) 

. Pressurized Water Reactor (PWR) Project Tech- 
nical Progress Report for the Period June 24, 1958 
to August 23, 1958, WAPD-MRP-75. (Unclassified 
AEC report.) 

. Power Reactor Technology, 1(2): 30 (February 
1958). (Available from the Superintendent of Docu- 





ments, U. S. Government Printing Office, Washing- 
ton 25, D. C.) 


. G. V. Alm and M. H. Binstock, Aluminum-electro- 


cladding Studies for Proposed OMR Fuel Elements, 
NAA-SR-2704, Oct. 1, 1958. (Unclassified AEC 
report.) 


. Metallurgy Division Quarterly Report for April, 


May, and June 1957, ANL-5790, July 1958. (Un- 
classified AEC report.) 
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Graphite 


High-purity graphite is used extensively in 
reactors because of its low thermal-neutron 
absorption cross section, its relatively high- 
moderating power, and its good high-tempera- 
ture properties. These properties, along with 
its relatively low cost and ease of machining, 
make it attractive in many applications as a 
combined moderator and structural material. 
A detailed review of the properties andirradia- 
tion behavior of reactor-grade graphite can be 
found in references 1 through 12. A brief re- 
view, containing some of the more important 
and recent data on graphite, follows. 


Thermal! Neutron Absorption 
Cross Section and Density 


Because of the importance of graphite as a 
moderator, a large effort has been made to re- 
move all impurities in the production of reactor- 
grade graphite and thereby reduce its absorption 
cross section to a minimum. Commercial puri- 
fied graphite is produced in the United States, 
having a density of ~1.68 g/cm® and an absorp- 
tion cross section of ~0.0037 barn. In a 1958 
Geneva conference paper,'* the French reported 
their best production graphite has a cross 
section of ~0.0038 and a density of 1.68 g/cm’. 

These results compare favorably with the 
value for the absorption cross section of C” of 
0.0033 + 0.00015 barn (2200 meters per second) 
established by Hennig"! and the cross section of 
highly purified graphite of 0.0035 + 0.0003 barn 
(2200 meters per second) reported by Muehl- 
hause. "4 


Density 


The theoretical density of graphite has been 
measured at 2.27 g/cm*. However, reactor- 
grade bulk densities of only ~1.68 g/cm’ have 
been produced in large quantities. A special 
“double-extrusion” technique reported by the 
French'® increased the density of the graphite 


to ~1.76 g/cm’. This value is still much lower 
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than the theoretical density, and voids amount- 
ing to 20 to 30 per cent of the volume are in- 
cluded in the graphite. These voids consist 
largely of interconnected micron-size holes 
which make the graphite permeable to gases 
and some liquids. 


Compatibility of Reactor 
Coolants with Graphite 


Although inert to most materials, graphite at 
high temperature is subject to oxidation. The 
large porosity presents a large surface area for 
oxidation, and consideration must be given to 
this possibility if oxidizing gases, such as air 
or CO,, are used as coolants in contact with 
graphite. Air will oxidize graphite at relatively 
low temperatures, and the prime emphasis has 
been placed gn the use of CO, and helium as 
gaseous coolants. Helium, being inert, can be 
used to very high temperatures before the im- 
purity content of it presents a problem. Carbon 
dioxide, which is being used in the Calder Hall 
type reactors, apparently must be used at tem- 
peratures below 1000°F to prevent serious 
graphite oxidation. A recent Geneva paper” 
indicated that the oxidation losses of graphite 
from the Calder Hall reactors will amount to 
0.1 per cent per year. This paper also re- 
viewed some of the studies made by the British 
on the interaction of moderator graphite andthe 
CO, coolant in a radiation field. Under in-pile 
conditions it was reported that CO, andgraphite 
reactions definitely indicated that irradiated CO, 
reacted more rapidly with graphite than non- 
irradiated CO,, the rate of reaction being nearly 
independent of temperature but varying almost 
linearly with radiation intensity. Three tech- 
niques were presented which can be used to 
inhibit these reactions: (1) Gas phase cracking 
of hydrocarbons to produce a resistant coating 
on the graphite; (2) Impregnation with phos- 
phorous compounds; and (3) Deposition of a SiC 
coating. 

Studies have also been conducted inthe United 
States on coating of graphite with NbC, TaC, and 
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ZrC to permit the use of graphite in reactors 
operating at high temperature in corrosive 
atmospheres.!*:!? These coatings are applied by 
a metallic vapor deposition technique in which 
the rate of coating is controlled by the diffusion 
of carbon through the carbide. Sound, adherent 
uniform coatings were produced which withstood 
repeated heating and cooling. Extensive thermal 
cycling tests were not conducted. 


General Properties 


The properties of unirradiated graphite have 
been adequately covered in the literature (ref- 
erences 1 through 7) and are not reviewedhere. 
However, further data on the unirradiated prop- 
erties of graphite at high temperatures 
(>2000°F)* and at very low temperatures (as 
low as 4°K)'® have also been reported in the 
references noted. 


Irradiation Behavior 


Of major interest to the reactor designer is 
the effect of irradiation on the properties of 
graphite. Much information of this type was 
covered in the 1955 Geneva conference for an 


Thermal conductivity: 
Electrical conductivity: 
Mechanical strength: 
Stored energy: 
Hardness: 

Dimensional changes: 


Markedly decreased 
Markedly decreased 
Markedly increased 
Increased 
Increased 

Normally expands 


Typical changes in some of the properties of 
graphite irradiated for various exposures are 
shown in Table VII-1. 

Above the 5000 Mwd/adjacent ton exposure 
at 30°C, saturation of irradiation damage effects 
has been reached, and the annealing of defects 
is almost as rapid as their rate of production. 

Irradiation behavior tests of graphite at 
higher temperatures has indicated that at tem- 
peratures in the 750 to 930°F range the ir- 
radiation growth usually associated with graphite 
can change to a shrinkage with increased ex- 
posure. A pilot of the test data showeda maximum 
voluine change of approximately + 0.04 per cent 
at an exposure of ~400 Mwd/adjacent ton fol- 
lowed by a volume decrease with higher burn-up 
which reached a volume change of approximately 
—0.2 per cent at 4000 Mwd/adjacent ton. At 
4000 Mwd/adjacent ton the rate of change of 


Table VII-1 PROPERTY CHANGES IN CSF GRAPHITE AT 30°C 


(Data taken reference 20) 





Exposure (Mwd/at)* 





Symbol Property Orientation 0 1000 2000 3000 4000 5000 6000 
Co Crystal layer spacing, A 6.71 7.16 7.48 7.65 7.73 7.78 7.78 
L. Apparent crystallite 

size in Cy direction, A 500 275 125 75 50 25 25 
K Thermal conductivity uaa 0.26 0.0082 0.0066 0.0063 0.0063 0.0062 
at 25°C, cal/cm-sec-°C Parallel 0.40 0.0104 0.0078 0.0073 0.0071 0.0070 
SE Total stored energy, cal/g 0 285 435 550 610 630 
Total stored energy 
remaining after a 3 hr 
anneal at 1000°C, cal/g 0 60 160 250 300 330 





* Exposures are given in terms of megawatt days per adjacent ton. This unit is defined as the amount of reactor 
radiation received by the sample during the time required for the ton of uranium adjacent to the sample to generate 
one megawatt-day of fission energy. One Mwd/at is equivalent to approximately 3.9 x 10!" nvt thermal or 6.5 x 10" 


nvt total. 


irradiation temperature of 30°C and at rela- 
tively low exposures.®»* In a paper presented 
at the 1958 Geneva conference further data were 


presented on the long term irradiation effects 
at 30°C and at higher temperatures.”° 

In general, the major effects of irradiation 
on graphite were: 


length parallel to the extrusion direction was 
~-—0.025 per cent AL/1000 Mwd/adjacent ton, 
and the rate of change of length perpendicular 
to the extrusion direction was ~—0.015 per cent 
AL/1000 Mwd/adjacent ton. 

This effect was also noted in some funda- 
mental studies on graphite which were made to 
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determine the relation between internal struc- 
ture and dimensional stability.2’:?? The results 
of the studies indicated that the net dimensional 
change obtained at a given irradiation tempera- 
ture is dependent on two competing effects: (1) 
crystal expansion and (2) shrinkage of inter- 
crystalline voids. In tests conducted at tem- 
peratures of 660 to 850°F, it was noted that the 
rate of crystal expansion was reduced under 
irradiation, and shrinkage of intercrystalline 
voids caused a net contraction of bulk graphite. 
In the irradiation experiments reported in 
reference 20, graphite pore size distribution 
showed a change from a peak at a pore radius 
of 200 A for zero exposure and 750 to 930°F to 
a peak at ~120 A for 4500 Mwd/adjacent ton, 
indicating a shrinkage of the intercrystalline 
voids. Preliminary data reported’ on the effect 
of irradiation on graphite operating at tem- 
peratures above 1000°F show: some increase 
in crystallite damage; very little further in- 
crease in the contraction rate at 1380°F; a 
greater contraction rate at ~1800°F; and a 
decrease in the stored energy. 

A radiation damage effect of considerable 
importance in low-temperature reactors is the 
energy storage effect. This phenomenon is due 
to the strain energy caused by distortion of the 
graphite crystals under irradiation. This energy 
will appear as heat when the irradiated graphite 
is heated above its irradiation temperature 
during annealing of the graphite. Typical an- 
nealing spectra curves, obtained by Woods et al.® 
by raising the graphite temperature at a con- 
stant rate and measuring the heat energy pro- 
duced, showed a thermal energy release with 
the peak release occurring at an annealing tem- 
perature of ~400°F. Thus, if low-temperature 
(<200°F) irradiated graphite is raised rapidly 
to this temperature, self heating can occur. The 
storage of energy is a function of the irradiation 
temperature, and the peak energy release at 
~400°F disappears with irradiation tempera- 
tures much above 250°F. 

The effect of the spontaneous energy release 
was presented in reference 8, where it was 
calculated that a specimen irradiated at 30°C 
to exposure of ~5000 Mwd/adjacent ton and 
heated to about 200 to 450°F, could achieve a 
spontaneous increase in temperature to ap- 
proximately 1800°F. 

The Windscale accident occurred during a 
scheduled heating of the reactor graphite for 
the release of stored energy. It should be noted 


that the safe release of stored energy can be 
accomplished with proper precautions. An ac- 
count of the energy release program carried out 
on the BEPO reactor on Mar. 14 and 15, 1958, 
is covered in detail in a 1958 Geneva conference 
paper.”° 

Other studies on the effects of irradiation on 
the properties of graphite are covered ina 
series of papers presented in BNL-489."! Of 
particular interest was a paper on the effect 
of fission fragment damage to graphite. Such 
damage can occur when the reactor fuel consists 
of a dispersion of UO, in graphite (Kernohan, 
reference 14). Tests conducted with graphite—5 
wt.% UO, specimens containing UO, particles of 
four different particle sizes indicated that at 
low irradiation temperatures and exposures 
(30 and 88°C at 0.1 per cent U atom burn-up) 
the specimen with the largest particle size 
(586 ».) suffered the least damage. The amount 
of damage appeared to be inversely propor- 
tional to the diameter of the UO, particle. Both 
conventional and fission fragment damage was 
noted in all specimens, but the fission fragment 
damage could not be annealed out. 


Zirconium and Zirconium Alloys 


Because of the low thermal neutron absorption 
cross section and relatively good mechanical 
properties of zirconium and its alloys, this ma- 
terial is of interest for various reactor design 
applications. One of the most important prop- 
erties which determine whether use can be 
made of zirconium or its alloys is their com- 
patibility with the coolant. Recent data on the 
corrosion of zirconium and zirconium alloys 
in water, steam, and gas environments are 
reviewed below. 


Compatibility in 600°F Water 


In general, zirconium and especially the 
Zircaloy-2 alloy have shown good corrosion 
resistance in 600°F water. Battelle Memorial 
Institute,*4 in a brief review of progress in 
zirconium technology reported corrosion rates 
for zirconium and Zircaloy-2 as shown in 
Table VII-2. 


Compatibility in 9OO°F Steam 


Zirconium and its alloys are known to have 
much poorer corrosion resistance in steam 
than in water at high temperatures. In an 
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Table VII-2 CORROSION RATES IN WATER 
AT 600°F 





Pre-breakaway 





Avg. weight Avg. 
gain, penetration rate, * 
mg/dm?/mo mils/yr 
Zirconium 4.2 0.088 
Zircaloy-2 5.4 0.113 


Post-breakawayT 


Zirconium 100 2.09 
Zircaloy-2 7.0 0.146 





* The average penetration rate is calculated on the 
basis of penetration of the original metal under uni- 
form corrosion conditions over the entire specimen 
surface due to conversion of zirconium to ZrQ,. 

tOccurs after many thousands of hours exposure 
to 600°F water. 


investigation of the problems involved in the 
development of a zirconium base alloy which 
would be corrosion-resistant in superheated 
steam, Nuclear Metals, Inc. (NMI) has made a 
number of experiments with zirconium and 
zirconium alloys exposed to 900°F, 1500 psi, 
superheated steam.”> The results of these ex- 
periments showed that zirconium and annealed 
samples of Zircaloy-2 and Zircaloy-3 undergo 
flaking of the oxide within one day of exposure 
to the 900°F steam. The tin content of the 
Zircaloy alloys was reported to be very det- 
rimental to the corrosion resistance to steam. 
Alloys of zirconium containing a minimum of 
0.2 wt.% nickel, or ~1 wt.% chromium or 
~2 to 3 wt.% iron were tested by NMI and 
found to have better corrosion resistance than 
Zircaloy-2; flaking in the 900°F steam en- 
vironment did not occur for approximately 40 
to 60 days. In all of the alloys tested, spalling 
or flaking occurred ata corrosion product weight 
gain of approximately 120 to 160 mg/dm’. These 
alloys still do not represent usable materials 
in a superheated steam atmosphere although 
they are considerably better than pure zir- 
conium, Zircaloy-2, or Zircaloy-3. 


Compatibility with a He-CO2 Gas Mixture 


In a study of the effect of reactor atmosphere 
on zirconium and zirconium alloys Hanford’® 
has reported on the in-reactor corrosion rates 
vf zirconium and Zircaloy-2 exposed to a 50 
vol.% helium-50 vol.% CO, atmosphere (Table 


Table VII-3 CORROSION OF ZIRCONIUM AND 
ZIRCALOY-2 IN A 50 VOL.% He — 50 VOL.% 


CO, GAS MIXTURE* 








Extrapolated Extrapolated 
weight gain, penetration, 
2 -3 ; 
Temp., mg/cm 10™ in. 
°C in3yr inlOyr in3yr inl0yr 

Zirconium 400 2.5 5.35 0.43 0.93 
Zircaloy-2 400 5.8 10.0 1.6 1.73 
Zirconium 500 9.4 20.0 1.56 3.5 
Zircaloy-2 500 32.5 89.0 5.5 15.3 





* Data taken from reference 26. 


Table VII-4 CORROSION OF ZIRCONIUM AND 
ZIRCALOY-2 IN DRY CARBON DIOXIDE AT 120 PSIA* 








Extrapolated Extrapolated 
weight gain, penetration, 
eis mg/cm? 1073 in. 

Material °C inlyr in3yr inlyr in3yr 
Iodide Zr 400 0.29 0.48 0.05 0.085 
C-melted Zr 400 0.33 0.56 0.06 0.097 
Arc-melted 

Zril 400 0.50 0.96 0.087 0.17 
Arc-melted 

Zr 2 400 0.41 0.66 0.071 0.11¢ 
C-melted 

Zircaloy-2 400 0.52 0.8 0.090 0.14 
Iodide Zr 500 11.0 31.0 1.91 5.4 
C-melted Zr 500 33.0 120.0 5.75 21.0 
Arc-melted 

Zri1 500 27.0 149.0 4.7 24.4 
Arc-melted 

Zr 2 500 23.0 100.0 4.0 17.4 
C-melted 

Zircaloy-2 500 130.0 450.0 22.6 78.0 





* Data taken from reference 27. 


VII-3) at a weighted average temperature of 
430°C for 886 days. The results of the tests 
are compared with British data?’ on the cor- 
rosion of zirconium and Zircaloy-2 in a dry 
CO, atmosphere in Table VII-4. 

Comparison of the two tables shows that a 
much lower corrosion rate at 500°C was ob- 
tained by the Hanford tests. Because of the dif- 
ferences in the test conditions, the reason for 
the lower corrosion rate is not apparent. Hanford 
indicates that the decrease was primarily due to 
a delay in the inception of breakaway corrosion 
caused either by the irradiation or some con- 
stituent in the gas mixture; impurities of N,, 
H,, CO, A, and O,, in the gas mixture were said 
to sometimes exceed 1 per cent. 
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Miscellaneous Components 


The Radiation Effects Information Center, 


Battelle Memorial Institute, recently released 
information on the effects of nuclear radiation 


on various electronic components, 
glass, and hydrocarbon fuels. 


gaskets, 
28-32 
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FUEL CYCLES 





u*® Problem in Isotope 
Separation Plants 


If the fuel of an enriched-uranium reactor is 
recycled a number of times through a diffusion 
plant for re-enrichment, the U6 content may 
build up to a rather high value because the rate 
of burnout in the reactor is low and there is no 
other effective sink for the U’** produced except 
the waste stream from the separation plant. The 
radioactivity and neutron absorption problems 
connected with the U?** (and the higher isotopes 
produced from it in the reactor) have been con- 
sidered in reference 1, which was discussed in 
a previous issue of this Review.’ Reference 3 
considers the effect of the U’** content on the 
cost of re-enrichment. The reference contains 
theoretical curves for a typical diffusion plant 
which show the equilibrium concentration of a 
in recycled fuel as a function of reactor feed 
enrichment, with the degree of burn-up per re- 
actor pass as a parameter. The U?** concentra- 
tion curves, for burn-up in the range 20 to 60 
per cent, show that for reactor feed enrichments 
in the range 30 to 50 per cent U?*> the U?** con- 
tent is comparable to the U®*> content, and that 
for higher enrichments U?** is the principal con- 
stituent other than U2*>, The cost increase con- 
sidered is only that due to the additional separa- 
tive work which is done on U**® in achieving the 
equilibrium concentration corresponding to the 
specified U**® content. This extra cost is of the 
order of 16 per cent for a reactor feed of 80 
per cent U5 content, decreases with U?*> con- 
tent, and is less than 1 per cent for U?** con- 
tents smaller than about 10 per cent. The 
equilibrium U?** content and the additional 
separation costs are greatly decreased if some 
of the diffusion plant output is bled off for other 
purposes such as fueling new reactors. A diver- 
sion of 25 per cent of the output decreases the 
u*** content of 50 per cent enrichment fuel by 
about a factor of 5. Itthus appears that the extra 
separation cost resulting from U** is not im- 
portant in low-enrichment reactors andis likely 
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to be small even for high-enrichment reactors 
if there is some diversion of the diffusion plant 
output. These conclusions apply only if the 
equilibrium U?** content, which corresponds to 
the diffusion plant output of fuel having the speci- 
fied U?*® content, is tolerable. The neutron ab- 
sorption and radiation problems connected with 
this equilibrium concentration still remain. 


Isotope Build-up in Recycled Fuel 


The economic operation of a reactor on the 
thorium cycle requires the recovery and use of 
the U?*? produced during the fuel irradiation. 
The buildup of extraneous isotopes in recycled 
fuel has been discussed, for both the uranium 
and thorium cycles in reference 4, which was 
reviewed previously.’ A recent report by St. John 
and Toops® discusses a number of factors re- 
lating to the accumulation of U’ in the U?* 
produced during irradiation of thorium andpre- 
sents a number of details in the derivation of 
mathematical expressions related to the forma- 
tion of U?8*, Uranium-232 is the primary source 
of gamma-ray activity in recycled U*** and, 
through its daughter Th’*, is a contributing 
factor to the increased gamma-ray activity of 
the salvaged thorium. The decay of @-active 
u*??_ (74-year half-life) leads to a number of 
radioactive daughter products, and it is the 
energetic gamma radiation from these daughter 
products which complicates the handling and 
fabrication of fuels using recycled U***.* Since 
the 1.9-year Th’”® is one of the daughter prod- 
ucts of U’**, the accumulation of this isotope, 
prior to chemical separation of the spent thorium 
fuel element, contaminates the thorium re- 
covered in the chemical processing stages. This 
contamination can be decreased by processing as 
soon after irradiation as possible, soas to mini- 
mize the growth of Th?”°, 

Uranium-232 is formed during the irradiation 
of thorium by either of two schemes: 





*ThC (Bi?!) emits a 2.2 Mev gamma ray. 
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Since the cross section for the n,2n processes re- 
quires high energy neutrons and would be ex- 
pected to be approximately the same for Th?” 
or U?*?, the rate controlling step of the processes 
above will be the n,y reaction on Th?” or Pa??!, 
The amounts of U’* formed by the two processes 
would then be in the ratio of n,y cross sections 
and would therefore be predominately (~97 
per cent) the second of the two schemes above. 
The relatively large capture cross section of 
u’*2 would tend to compensate for the U*®” 
production in the first scheme, i.e., the n,2n 
reaction on U”*?, 


As sources of fast neutrons for the n,2n reac- 
tion reference 5 considers both the neutrons 
originating in the thorium lump itself (from U2** 
fission) and those originating in fuel elements ata 
distance. The computation of the U?*? production 
by these neutrons is set up for solution by the 
Monte Carlo technique on anIBM-650 computer. 
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HEAVY-WATER REACTORS: LEAKAGE DETECTION 





The D,O moderated reactor is under develop- 
ment in a number of countries; interest in the 
reactor type is prompted mainly by the high 
neutron economy attainable and by the possi- 
bility of achieving rather long fuel life withura- 
nium of natural enrichment (see pp. 7 to 22.) 
Power reactors cooled and moderated by D,O 
are under construction or active development in 
Canada,' Sweden,” Norway,’ and the United 
States.‘*® D,O-moderated reactors cooled by 
gas’ and by sodium® are under development in 
the United States, and a large gas-cooled D,O- 
moderated reactor is reported to be under con- 
struction in Czechoslovakia.‘ 

Perhaps the most important disadvantage of 
D,O is its high cost, which in the United States 
amounts to $28 per pound. Obviously pertinent 
questions in the evaluation of D,O reactors are 
the probable rate of loss of D,O, the rate of 
contamination of D,O by H,O, and the costof re- 
enriching contaminated D,O. Babcock‘ has stated 
that over a two-year period the D,O losses from 
the Savannah River production reactors have 
averaged about 3 per cent per year. Hehas also 
estimated that an average entry rate of about 
2 lb per day of H,O into the D,O might be ex- 
pected in a large power reactor. As a means of 
estimating the cost of removing the leakage H,O 
from the D,O, he cites the information that 1000 
lb per hour of steam used to heat the calandria 
of a well-designed vacuum distilling tower will 
remove 1 lb of H,Oper day froma 99.75 per cent 
D,O liquid asa 75 per cent D-25per cent H mix- 
ture. If heat from the steam is assumed to cost 
30 cents per million Btu, then the heat cost for 
removing H,O would be about $8.20 per pound. 
Each pound of removed H,O would be accom- 
panied by 3 lb of D,O; the 75 per cent D,O-H,O 
mixture would be sent to a central D,O plant 
for re-enrichment at a cost of perhaps $3to $4 
per pound of contained D,O, or $9 to $12 
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per pound of H,O. Thus the operating cost for 
H,O removal would amount to about $20 per 
pound of H,O. To this must be added capital 
charges against the cost of the distilling tower. 
On this basis, it would appear that for reactors 
equipped with their own distilling towers the 
cost of H,O inleakage, per pound of leakage, is 
comparable to the cost of D,O outleakage. 

The foregoing estimates, of course, were made 
for reactors cooled and moderated by D,O. The 
leakage rates which characterize specific plants 
will depend entirely on the plant designs. A 
number of reactor designers hope to make con- 
siderable improvements over the estimates. 
For reactors which use D,O only as moderator, 
the leakage expectation should be much lower, 
since the moderator may be unpressurized and 
does not require rapid circulation, andaccess to 
the moderator is not generally required for fuel 
changes. 

Regardless of anticipated leakage rates it is 
obvious that effective monitoring for both in- 
leakage and outleakage is demanded by economic 
considerations. 

Another area of concern related to D,O leak- 
age is the protection of personnel. D,O that has 
been irradiated in a high neutron flux is diluted 
to a predictable extent by tritiated water (T,O or 
DTO), tritium being formed from deuterium by 
neutron absorption. Listed below are the three 
main nuclear reactions related to tritium: 


jH? + gn! — ,H®? + y [o = 5.7 x 10 barns] 

iH? _12-4 "| jHe® +_,e° [18 kev] 

»He® + on' — ,H'+ ,H® [o = 5400 barns] 
Leakage of tritium gas or tritiated water 


vapor to the air in a personnel access area pre- 
sents a significant health hazard because of the 
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long half-life of the beta emitting decay of 
tritium. Tritiated water vapor is readily ab- 
sorbed by lung tissue or skin and becomes uni- 
formly distributed throughout the body. The 
following are significant data concerned with 
the hazard: 


1.3 mc (ref. 8) 
10 mc (ref. 9) 


Maximum permissible 
body burden: 

Maximum permissible 
concentration in 
atmosphere: 

Maximum permissible 
concentration in 
water: 

Effective or biological 
half life in body: 


5 x 10 *uc/cc (ref. 8) 
2 x 107° ue/cc (ref. 9) 


0.2 uc/cc (ref. 9) 


19 days (ref. 9) 


Many types of serious leakage to the air may 
go unnoticed with inadequate instrumentation. 
These include: vapor leakage from shaft seals 
and gaskets; leakage through small flaws ina 
hot pressurized system; helium leakage (many 
reactor designs employ a helium blanket to pres- 
surize the D,O and flush off fissiongases or de- 
composed water); D,O removal with fuel ele- 
ments and transferral to handling equipment 
surfaces and to cooling canal and retention basin 
water (and also to the air above such facilities). 

In the following paragraphs various instru- 
mentation approaches to the separate phases of 
the leak monitoring problem are reviewed. No 
attempt is made to cover all possible types or 
variations or to propose an integrated system. 


Quantitative Analysis 
of Water Mixtures 


Infrared spectrometry isa practical andaccu- 
rate technique for quantitatively analyzing mix- 
tures of liquids. This technique is applicable 
where the components of the mixture have 
recognizably different absorption characteris- 
tics for infrared radiation. 

A study of D,O leak-detection methods basedon 
infrared spectrometry, by Stevens and Bayly” at 
Chalk River resulted in the development of 
methods for determining accurately the propor- 
tions of D,O, H,O, and HDO in mixed water. 
These liquids fortunately have characteristic and 
differing infrared radiation absorption proper- 
ties. The work was primarily concerned with the 
very high andvery low D,O concentration ranges. 


D,O strongly absorbs infrared radiation of 
2.94 » wave length. HDO and H,O do not absorb 
this wave length as readily. Therefore, the de- 
gree of transparency of a mixture to this wave 
length can be taken as an indication of the D,O 
concentration. Using a single beam spectrome- 
ter, properly calibrated with known mixtures, it 
is possible to analyze a water mixture and de- 
rive the per cent concentration of D,O witha fair 
degree of accuracy. However, by splitting the 
light source into two equal intensity beams of 
the same wave length, passing one through the 
sample as above while the other passes through 
a known reference absorber, the transmission 
or absorption difference may be measured. 
Variations in the sample beam then result in 
magnified variations inthe difference indication. 
A commercial double beam spectrometer has 
been adapted for use in this manner. A soft 
glass plate is used as the reference absorber 
in this arrangement. In the range of D,O con- 
centrations from 99.5 to 99.9 per cent, and with 
proper control, precision of + 0.003 per cent D,O 
is reported. Analysis time is approximately 
7 min per sample. 


A difference technique similar to that de- 
scribed above is used to measure the degree of 
dilution of H,O by HDO and D,O." In this case 
the width of the H,O absorption band centeredat 
3.1 H would obscure the D,O band at 2.94 y. 
However, HDO has a usable absorption band at 
4.0u, and, as HDO is normally the major diluent 
in this instance, it can be successfully used. A 
cell filled with natural distilled H,O is used in 
the reference beam. After preparation ofa cali- 
bration curve from known samples, precision of 
+0.01 per cent HDOor D,Ointhe 1 per cent con- 
centration range, to +0.003 per cent in the 
0.03 per cent range, is expected. Control of cell 
temperature in the sampling and calibration pro- 
cedures is critically necessary for accuracy. 


Continuous -flow 





Heavy-water Analyses 


An instrumentation system designed for con- 
tinuous monitoring of D,O concentration in the 
Chalk River establishment is also described 
in reference 10. This system makes use of a 
double beam nondispersive analyzer developed 
by the Perkin-Elmer Corp., under the trade 
name of “Tri-Non.” This system uses, as its 
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primary sensing device, what may be calleaa 
“microphone-condenser” type of transducer. 
This device consists of a double-walled gas 
chamber in which the inner wall is flexible. The 
inner and outer volumes of the chamber are 
joined by an orifice so that rapid changes in 
pressure in the inner chamber will result inin- 
ward and outward flexure of the inner walls with 
the orifice maintaining a long-term average 
pressure balance between the volumes and also 
acting as a pneumatic filter to prevent oscilla- 
tion. One plate of an electrical condenser is se- 
cured to the flexible wall, while the other plate 
holds a fixed position in the outer volume of 
the chamber. Nondispersed infrared radiation 
passed through transparent end windows in the 
gas chamber is absorbed by the gas, causing 
heating and expansion in proportion to the 
absorptivity of the selected gas for the particu- 
lar wave length of radiation used. In the sys- 
tem discussed in the reference, a mechanical 
“chopping” device is used to alternately permit 
the radiation from a beam passing through a 
variable optical attenuator, and a beam passing 
through a sample cell, to enter the gas chamber. 
Preliminary matching of the attenuator absorp- 
tion to the absorption of a D,O sample of known 
concentration, inthe sample cell, results in equal 
absorption by the gas cell from each beam as 
they are alternated. Therefore, the condenser 
plates remain at constant spacing and the elec- 
trical capacitance is unchanged. Now, by passing 
a continuous unknown sample stream through the 
sample cell and continuing to alternate the beams 
at a 6%, counts/sec rate, the heating effect and 
the capacitance will vary in proportion to the 
difference between the beam intensities at the 
selected wave length of radiation. The varying 
capacitance is used inanelectrical circuit which 
develops an a-c signal proportional inamplitude 
to the absorption difference. This signal is then 
used in a null balance servo system which 
operates to continually match the variations in 
sample transparency by adjusting the variable 
optical attenuator and also to develop an output 
signal proportional to the per cent of variation 
from the reference condition. This latter signal 
may be continuously recorded on a strip chart. 
Analyzing for absorption in the 3 band, for 
changes in D,O concentration over the 99.5 


to 99.9 per cent range, a change of 0.002 per 
centis said tobe readily detected. The response 
time is of the order of 3 sec. Unfortunately, the 


analyzer is very sensitive to sample tempera- 
ture. Therefore, provision must be made to 
control the sample stream temperature to with- 
in + 0.5°C of the calibration temperature before 
entry into the analyzer. 


Heavy-water Vapor 
Leakage Detection 


A semiportable device for probing for the 
point of D,O leakage, and for assessing the leak 
magnitude for D,O vapor leakage to the air, has 
been devised by Stevens and Bayly and is re- 
ported in reference 10. This device also makes 
use of a Tri-Non analyzer as its basic sensing 
element. The specification for the design of the 
system includes the requirement for minimum 
sensitivity to varying quantities of H,O vapor, 
CO,, and typical hydrocarbon vapors. Successful 
operation is reported using D,O vapor in the 
sensitive gas chamber, and a mixture of the 
gases to be discriminated against in the filter 
assembly of the reference beam attenuator cir- 
cuit. The air to be analyzed is drawn into the 
instrument at a rate of 75 cm*/sec through a 
pistol gripped probe, dust filter, and flexible 
hose. A meter is provided on the probe, which 
is calibrated in grams D,O per day at sampling 
rate. The device is also equipped to generate an 
audible signal whose pitch varies with D,O con- 
centration. Response time is approximately 5 
sec, and the instrument is reportedto be capable 
of sensing from 1 part D,O vapor in 6000 parts 
air to 1 part in 30 parts (approximate saturation 
at 25°C). 


Monitoring for Tritium in the Air 


Care must be taken that the presence of 
tritium gas or tritiated water vapor (T,0, HTO) 
in the air is noted rapidly. Therefore, in in- 
stallations where heavy water is subject to neu- 
tron irradiation, a tritium monitoring system in 
personnel areas is provided. 


As the characteristic decay beta emission has 
low average energy (18 kev), ordinary window 
type beta counters are not generally suitable, 
due to the extreme thinness required of a window 
for this energy. Most devices for tritium moni- 
toring therefore use ion chambers into which the 
suspected sample is directly or indirectly 
inserted. 
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One such device is briefly described in refer- 
ence 11. This small, portable device, designed 
for simplicity and rapid response, provides for 
drawing in a continuous flow of air through an 
electrostatic chamber, where the ion current, 
due solely to the ionization of the air by the 
weak beta emission, is sensed by anelectrome- 
ter tube. Current amplification is provided to 
yield an output sufficient to operate a sensitrol 
relay, which in turn sets off an alarm when full 
scale is reached on any one of the four ranges 
of the instrument, covering full-scale ionization 
currents from 107" to 10~* amp. This device is 
somewhat ambiguous in operation, as it senses 
ionization caused by other than tritium sources. 
However, its rapid response to a change in 
tritium background level warrants its investi- 
gation for general area monitoring. 


A device of interest’ developed at the Bureau 
of Standards for spot sampling of air for 
tritiated vapor is stated to be accurate for 
measurement of HTO vapor concentrations as 
low as that providing the maximum permissible 
exposure on a continuing 24-hr basis (e.g., 2 x 
10-° yc/ml).2 The device uses an ionization 
chamber of 250 cm® capacity, and a vibrating 
reed electrometer for current measurement. 
In use, the chamber is evacuated in the labora- 
tory, carried to the area to be examined, then 
opened to the air. The discharge of the elec- 
trometer is timed and referred to the calibra- 
tion curve from which the HTO concentration 
may be read directly. 


Routine monitoring for tritium in the plant 
atmosphere, as a background control means, is 
also accomplished by collection of atmospheric 
moisture from dehumidifiers or desiccators.® 
Standard ion chamber methods are used in the 
laboratory to evaluate the tritium content. 


In a more recently developed method, " triti- 
ated water is mixed with a liquid phosphor solu- 
tion. Tritium content is then measured by scintil- 
lation counting. 


A semiportable continuous air sampling de- 
vice which has been used extensively at the 
Savannah River Plant® for probing for a point of 
leakage is the Kanne Chamber." This device 
consists of a large ion chamber, a pump to 
sweep air through the chamber, and the associ- 
ated electronics for measuring the ion current, 
which after calibration indicates the tritium 


concentration. 


Monitoring for Tritium 
Contamination of Surfaces 


Tritium contamination of a surface is difficult 
to detect. The usual method of detection consists 
of wiping the suspected surface with filter paper, 
then inserting the contaminated paper into a 
counting chamber. It is understandable that this 
method cannot be considered quantitative under 
any conditions, because of surface roughness, 
porosity, and crevices that cannot be touched by 
the filter paper. Therefore, if a filter wipe 
(smear test) reveals some degree of contamina- 
tion, and the part cannot be removed for chemi- 
cal cleaning, it is not considered unusual to paint 
over the surface and thereby shield against the 
possibility of undetectable contamination in the 
area of an arrested D,O leak. 


However, Colvin'® at Savannah River Labora- 
tory has recently devised an open-ended ioniza- 
tion chamber, reported to be capable of quanti- 
tative measurement of tritium on rough or 
smooth surfaces removed as far as 7 in. from 
the open end. The chamber is a 3',-in.-diameter 
dished vessel, 1 in. deep, constructed of lucite 
and painted with conductive paint. A center elec- 
trode, insulated from the body, extends to within 
'’, in. of the plane of the open end. The chamber 
is operated as a dosimeter to integrate the ioni- 
zation current. It is charged to a specified volt- 
age and then placed upon the suspected surface. 
After a specified time has elapsed, the chamber 
voltage is remeasured. The decrease in voltage 
is proportional to the amount oftritium present. 
The chamber is charged and read by use ofa 
standard electrometer. The configuration of the 
electrostatic field between the center electrode 
and the shallow outer electrode, as stated above, 
extends vA in. beyond the end of the chamber; 
therefore the weak tritium betas, which have a 
range of only ie in. instandardatmosphere, can 
contribute to chamber ionization even though they 
emanate from crevices or other unevenness 
having a considerable depth. Using this device it 
is stated to be possible to measure as low as 
12,000 decays per minute, with a signal-to- 
noise ratio of 10. 
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BOILING REACTORS: TURBINE CORROSION AND EROSION 








An investigation of erosion and corrosion by 
wet oxygenated steam, by Cataldi, Cheng, and 
Musick, has been reported recently in the 
Transactions of the American Society of Me- 
chanical Engineering.’ With specific application 
to Dresden?** dual-cycle steam conditions, the 
article discusses a series of preliminary labora- 
tory tests evaluating the resistance of various 
materials for turbines and related components 
to erosion (impact, washing, wire-drawing) and 
corrosion (weight loss, pitting depths, crevice 
and galvanic attack) inthe wet, oxygenated steam. 
Materials tested included carbon and low-alloy 
steels and cast iron, austenitic and ferritic 
stainless steels, monel, cupronickels, leaded 
bronze, stellite, and Ductile Ni-Resist No. 3. 


Steam from a direct-cycle or dual-cycle 
boiling water reactor differs from usual power 
plant steam in several respects: it is radio- 
active, it may be wet, and it contains some free 
hydrogen and oxygen. These characteristics are 
discussed in references 4 and 5, which were re- 
viewed previously.® The General Electric Com- 
pany, in order to gain experience with the effect 
of oxygen on corrosion of metals in high- 
temperature water-steam mixtures, undertook 
an experimental program with the twofold pur- 
pose of: (1) aid in the selection of materials 
for the turbine of the Dresden Nuclear Power 
Station; (2) fundamental study of the effects of 
temperature, pressure, oxygen, velocity, and 
pH on erosion and corrosion by water-steam 
mixtures. Dresden cycle conditions are: (1) wet 
reactor steam at 1015 psia (545°F); (2) turbine 
inlet primary steam at 869 psia (525°F) and 0.5 
per cent moisture; (3) secondary steam at 415 
psia (450°F). 

The tests conducted are grouped under three 
headings, corresponding to the facilities used. 


Test 1. Erosion and corrosion with wet Sche- 
nectady Plant steam: 350 psia and 429° maxi- 
mum; pH 6.8; with 0.05, 80, and 110 ppm oxygen. 

Test 2. Erosion and corrosion with wet steam 
using desuperheated steam from 14th stage of a 
turbine at Mystic Station of Boston Edison Com- 


pany; 43 psia and 272°F, 70 to 110 ppm O,, pH 
8.8 to 9.5. 

Test 3. Pitting, crevice, and galvanic corro- 
sion tests inwater and steam inautoclaves: 1015 
psia (545°F), 50 and 100 ppm O, (initially), pH 7. 


Some of the results of the tests are given in 
Table X-1. For complete information onthe test 
specimens and conditions the original article 
should be consulted. 

The following summary and conclusions were 
given by the authors: 


1. On the basis of this preliminary evaluation, 
erosion and corrosion of aturbine in a boiling water 
nuclear power system can be maintained within 
tolerable limits with proper design and the selec- 
tion of special materials for critical areas. This 
evaluation is being supplemented by continuing labo- 
ratory tests,and by the operation of a 5000 Kw tur- 
bine directly connected to 210,000 Kwboiling water 
reactor at the Vallecitos Atomic Laboratory (Atomic 
Power Equipment Department). 

2. In wet steam with low (0.05 ppm) and. high 
(100 ppm) oxygen, stainless steels, Monel: and 
Stellite 6 were most resistant to ‘‘wire drawing’’ 
erosion. Ductile Ni-Resist 3, low nickel steels, and 
low chromium steels were better than low molyb- 
denum and carbon steels; plain cast iron was the 
least resistant material tested. 

3. In wet steam at pH 7 and 430°F, there was no 
effect of oxygen on erosion in tests with high and 
low oxygen concentrations, except for BTH leaded 
bronze. With low oxygen, BTH leaded bronze re- 
sembled low chromium steel in erosion; with high 
oxygen, it resembled plain cast iron. 

4. In wet steam at 300-430°F, and pH 7, with low 
and high oxygen, stainless steel, Stellite No. 1, and 
Ductile Ni-Resist No. 3were most resistant to cor- 
rosion. Monel and low-alloy steels were superior 
to carbon steel. BTH leaded bronze was the least 
resistant material tested. At lower temperatures 
(300-400°F) and low oxygen, the cupro nickels were 
comparable to the most resistant materials; at-high 
oxygen, they very rapidly approached the extreme 
corrosion rate of BTH leaded bronze. 

5. At 300-400°F, pH 7, increasing. the oxygen 
level from low to high increased the corrosion rate 
of all materials tested except stainless steel. Inthe 
case of carbon steel, this effect of oxygen on corro- 
sion diminished rapidly with increasing tempera- 
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Table X-1 CORROSION RATES IN WET STEAM 
(Mils per Year) 











Test 1 Test 2 Test 3 

Loop or run 1 2 1 2 1 2 A B 
Temperature, °F 430 430 400 400 300 300 272 272 545 545 545 545 
Pressure, psia 350 350 245 245 67 67 43 43 1015 1015 1015 1015 
Moisture, % 9 9 7 7 5 5 6 6 Sat. Sat. Sat. Sat. 
Velocity, ft/sec 30 30 30 30 45 45 90 90 Rocking Autoclave 
pH 6.8 6.8 6.8 6.8 6.8 6.8 9.5 8.8 7 i 7 7 
Oxygen, ppm 0.05 110 0.05 110 0.05 110 70 110 50 100 ©6550 100 
Duration 1000 hr 1150 hr 1000 hr 1150hr 1000hr 1150hr 780hr 1150 hr 6 wk 6 wk 10 wk 10 wk 
Carbon steel: 

ASTM A106, A212, A201 2.0 2.2 1.0 1.4 0.2 1.0 2.6 16 5] SS tt 

ASTM A27 cast 2.5 CSOT Be Be 
Low-alloy steel: 

2¥, Ni, y, Mo, forged 0.9 1.3 0.4 0.6 0.3 0.4 a7 e- RF’ 88 ‘G4 

2, Cr, 1 Mo, (A155) 1.1 1.4 0.5 ee 0.2 0.4 0.6 0.8* 0.8* 0.5* 0.3* 
Cast iron: 

ASTM A48 1.9 22:30 @& (82 

Ductile Ni-Resist No.3 0.2 0.4 0.2 0.2 <0.1 <0.1 0.7 0.1 0.1 
Stainless steel: 

ASTM type 304L % Sen. 0.2 0.1 0.1 0.1 0.1 <0.1 0.9 0.2 9.2 

ASTM type 347 <0.1 0.1 Nil Nil Nil Nil Nil Oe O62 8 AA 

ASTM type 403 0.2 0.2 0.2 0.1 0.1 <0.1 Se 08 CE “O72 
Nonferrous alloys: 

90/10 Cu-Ni 0.4 4.6 0.2 0.6 1.2 

80/20 Cu-Ni 0.3 4.4 0.1 0.6 

70/30 Cu-Ni 0.2 6.8 0.1 1.3 

Monel 0.8 0.9 0.2 0.4 0.1 <0.1 0.7 

BTH Pb bronze 3.2 14.5 2.4 4.0 3.3 

Stellite No. 1 <0.1 0.1 <0.1 <0.1 <0.1 <0.1 





*2 Cr cast steel. 


ture. With the cupronickels, the effect of oxygen 
was much more severe at 400°F than at 300°F. 

6. Corrosion rates of all the materials tested, 
at 300-400°F, pH 7, high and low oxygen, increased 
with temperature, the logarithm of the corrosion 
rates varying approximately linearly with the re- 
ciprocal of the absolute temperature. 

7. At 430°F, pH 7, with high oxygen, there was 
no intergranular corrosion of sensitized types 304 
and 306L stainless steel specimens in 1000 hours, 
nor were there any cracks around stamped crosses. 

8. With low oxygen, at 300-430°F, pH 7, there 
was no pitting observed except with carbon steel, 
where the effect was minor. With high oxygen, all 
materials pitted especially at crevices. Stainless 
steel and Stellite No.1 were most resistant. Monel, 
cupronickel, and BTH leaded bronze were better 
than low-alloy steel. Carbon steel pitted most 
severely of the alloys tested. 


9. In autoclave tests in steam-water-oxygen 


mixtures (545°F, 1015 psia, 50 and 100 ppm O,) run 
for 6 and 10 weeks, the rates of corrosion decreased 
with time. The maximum depth of pitting, however, 
did not increase significantly. 


a 


10. Inthe autoclave tests, the order of resistance, 
as measured by the maximum depth of pitting, was 
the same as forthe resistance to general corrosion. 
There was less difference in the maximum pitting 
depths among the ferrous alloys tested than there 
was in the corresponding weight losses. This re- 
flects the greater number of pits in the least re- 
sistant alloys. 

11. There were no significant galvanic or crevice 
effects inthe autoclave tests. Naturally-formed de- 
posits of iron oxides were more effective in pro- 
moting pitting on all classes of steel tested than 
were artificial 0.002-inch crevices. 

12. Cast iron and low-carbon cast steel were 
subject to occasional very deep pits in the auto- 
clave tests. These may be related to inclusions 
and/or pin-hole porosity. 
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